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Ensaios de fundações rochosas de barragens 


e seu comportamento 


SUMÁRIO 


Este trabalho apresenta os resultados dos esforços 
de vários anos para interpretar ensaios de deformadbili- 
dade de fundações rochosas de barragens de betão. Os 
resultados de ensaios de deformabilidade «in situ» de- 
pendem da direcção e profundidade dos pontos ensaiados 
e das condições geológicas do local, sendo muito claras 
as correlações com a profundidade e inclinação das for- 
mações. São discutidos os métodos de interpretação 
usados e feitas comparações com resultados obtidos em 
laboratório, sendo apresentados os resultados de diver- 
sos ensaios «in situy em diferentes rochas. Indicam-se 
novos factos sobre os efeitos da injecção na deformi- 
dade agora encontrados e dão-se informações recentes 
sobre o comportamento das fundações em rocha, basea- 
das em observações com medidores de extensões de 
grande comprimento (extensometros). 


INTRODUÇÃO 


No projecto de grandes barragens de betão, é de 
importância básica o conhecimento da deformação da 
fundação devida às cargas impostas pela barragem, pela 
pressão da água intersticial e por outros efeitos. Um 
conhecimento da resistência ao corte da fundação é tam- 
bém essencial, especialmente ao longo das superfícies 
de resistência minima. Tem havido progressos nesses 
dois campos experimentais da Mecânica das Rochas, 
mas os métodos de ensaio mais seguros continuam a ser 
semelhantes âqueles que se desenvolveram há duas dé- 
cadas. Como exemplos temos os ensaios com macacos 
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SUMMARY 


The paper presents results of many year's efforts 
to interpret deformability tests of foundation rocks of 
concrete dams. Results of «in situ» deformability tests 
depend on the direction, depth of tests and the geologi- 
cal details of the site. Correlations with depth and dip 
and strike of the formations are very clear. Methods of 
interpretation used and comparisons with laboratory re- 
sults are discussed and results of many tests at site in 
various rocks are presented. New findings on the effects 
oi grouting in deformability are indicated and the most 
recent information is given on rocks foundation heha- 
viour, based on cbservation with long length strain me- 
ters (extensometers |). 


dentro de galerias e os ensaios de compressão triaxial 
em laboratório. Contudo, tais métodos tornaram possivel 
a descoberta de muitos aspectos do comportamento me- 
cânico de maciços rochosos, que muito contribuiram 
para os conhecimentos neste campo. No caso das deter- 
minações de deformabilidade, em granitos, os primeiros 
resultados indicaram variações muito importantes com o 
grau de alteração (que no geral diminui com a profun- 
didade) e com o grau de compressão natural a mostra- 
ram também variações com a direcção, principalmente 
em ensaios no laboratório de amostras do material ro- 
cha ('). Tornou-se então evidente que era necessária uma 
interpretação mais pormenorizada dos resultados. 
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Os ensaios tradicionais «in situr com macacos, em 
galerias, usados pelos autores para determinar a defor- 
mabilidade, incluem ciclos de carga sucessivamente 
crescentes, até uma tensão normal máxima superior à 
imposta pela barragem às suas fundações e a repetição 
dos ciclos de carga e descarga até esta tensão, Geral- 
mente, o «módulo de elasticidade» do maciço rochoso, 
«E», é tomado como a relação entre essa tensão máxima, 
oc e a média das deformações máximas, à, observadas 
nos ciclos repetidos, durante a carga, no centro da área 
carregada ("): «E» = od (1 — »') / 8, onde d é o diâmetro 
da área em carga e » o coeficiente Poisson da rocha. Este 
«módulo de elasticidade» é calculado supondo que o 
corpo é um sólido semi-infinito, elástico, continuo, homo- 
géneo e isotrópico. A razoável aproximação de tais con- 
siderações tem sido verificada, mas pode ser posta em 
dúvida ('), ('). Considera-se que a justificação para estas 
afirmações terá de vir também da observação directa do 
comportamento das fundações das barragens. 


Todos estes anos de estudo da Mecânica das Ro- 
chas têm conduzido a uma visão coerente do comporta- 
mento dos maciços rochosos. Ao estudar os diagramas 
de tensão-deformação em ciclos de carga e descarga 
dos ensaios «in situ» e ao considerar a existência de fis- 
suras finas, foi-se conduzido a propôr para os maciços 
rochosos um «modelo de esponja» que explica: (a) o 
assentamento progressivo observado nos ciclos suces- 
sivos de carga e descarga: (b) o aumento do «módulo 
de elasticidade» tangente com carga crescente até uma 
determinada tensão (quando começa a ruptura); (c) o 
«módulo de elasticidade» tangente muito mais alto 
quando começa a descarga, etc. 

Os resultados dos ensaios triaxiais em laboratório (') 
têm também posto em relevo não só os aspectos (b) e 
(tc), mas ainda que a deformabilidade de uma amostra 
para a mesma tensão axial diminui bastante com a pres- 
são confinante, o que é também considerado um efeito 
do fecho das microfissuras. Em tais ensaios, o «módulo 
de elasticidade», «E», é considerado como a relação 
entre a tensão axial total e a extensão longitudinal cor- 
respondente, apesar de se saber que se trata de sólidos 
anisotrópicos. 

Em comunicações recentes (*), ('), ('), os autores 
têm indicado os esforços feitos no sentido des uma in- 
terpretação racional dos resultados dos ensaios «in situ» 
e de encontrar maneiras mais eficientes para determinar 
a deformabilidade dos maciços rochosos. O propósito 
principal deste trabalho é o de relatar a continuação 
destes esforços, apresentar resultados de ensaios recen- 
tes e indicar valores provenientes da interpretação de 
medições de extensões em fundações das barragens 
em serviço. Também se apresenta um relato pormenori- 
zado do método de interpretação dos ensaios, agora 
usado. 

Entre os resultados mais importantes obtidos estão 
os que indicam, sistematicamente, para todos os locais 
de barragens, o aumento dos «módulos de elasticidade» 
dos maciços rochosos com a profundidade. Torna-se 
claro que a descompressão, a alteração da rocha e a 
fissuração mais intensa junto da superfície são factores 
primários determinando valores baixos dos «módulos de 
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elasticidade» dos maciços rochosos nesta zona. Contudo, 
a deformabilidade de um dado maciço rochoso num de- 
terminado local parece atingir um limite a partir de uma 
certa profundidade, se não houver modificação no grau 
do alteração, embora não seja igual para todas as direc- 
ções. À composição mineral, a matéria pétrea, a poro- 
sidade e as microfissuras da rocha também influenciam, 
evidentemente, as propriedades do maciço rochoso. 
Mesmo em granitos se tem notado que os planos das 
microfissuras têm uma orientação preferencial no mate- 
rial rocha ('). 


Sendo os maciços rochosos meios discontinuos 
(devido à existência de fissuras, falhas, xistosidade, 
planos de estratificação, etc.), como resultado imediato, 
o material tem um comportamento anisotrópico especial, 
As condições de sedimentação, as tensões tectónicas 
desviatórias e variáveis, o metamorfismo, a fracturação 
devida a tensões internas, sismos ou explosões, a recris- 
talização e outros factores também contribuem, de modo 
muito importante, para a variação das propriedades com 
a direcção, Os autores têm frequentemente chamado a 
atenção para o facto dessa anisotropia não poder ser 
assimilada a qualquer dos modelos matemáticos melhor 
estudados de materiais anisotrópicos da Teoria da Elas- 
ticidade. As tentativas de interpretação dos dados ob- 
tidos, quer em ensaios de campo em galerias quer no 
laboratório, supondo, por exemplo, um comportamento 
ortotrópico, monotrópico, ou isotrópico transverso, têm 
falhado nos casos de formações xistosas ou de cama- 
das sedimentares ('), ('), Os resultados das medições 
de extensões em diversas direcções em várias amostras 
tiradas em ângulos diferentes de um local e ensaiadas 
triaxialmente não estão de acordo uns com os outros 
(são incompatíveis se a rocha se suposer continua), in- 
dicando que as caracteristicas anisotrópicas duma mesma 
rocha mudam conforme a amostra, a direcção e o tipo 
das cargas ('). Este problema é esclarecido neste tra- 
balho. 


Em tentativas de interpretação da anisotropia dos 
maciços rochosos a partir dos resultados de ensaios «in 
situr na mesma zona, os cálculos baseados na isotropia 
transversal levaram à conclusão que os valores dos mó- 
dulos de elasticidade que poderiam ser deduzidos são 
muito próximos dos calculados ao assumir isotronia per- 
feita. Isto deu uma razão para aceitar os cálculos tradi- 
cionais dos resultados dos ensaios «in situp atrás refe- 
ridos ('). 

Considerando os resultados dos ensaios menciona- 
dos, parece muito importante fazer notar que as constan- 
tes de deformabilidade e os parâmetros de anisotropia 
dum maciço rochoso não dependerão apenas do próprio 
maciço (material, textura, estado de tensão, grau de des- 
compressão devido à excavação, etc.), mas também das 
tensões impostas ao maciço durante os processos de 
carga e descarga e a sua sucessão. 

Para as diversas direcções da carga, a possibili- 
dade de deslizamentos, ao longo de superfícies de mi- 
crofissuras, planos de xistosidade, ou planos de cliva- 
gem de cristais orientados, determina que a deformabi- 
lidade máxima ocorra quando a direcção da carga e o 
ângulo de atrito interno dessas superfícies permitem 
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cesso deslizamento. À influência mecânica de uma falha 
ou de uma família de diaclases pode também ser, em 
grando escala, semelhante aos efeitos das microfissuras, 
Nesta trabalho mostra-se que a orientação das microfis- 
suras e, esjecialmente, o deslizamento que nelas pode 
ocorrer, é de importância primária no comportamento 
«anisotrópicor observado em amostras de laboratório. 
À geometria do corpo carregado é evidentemente básico 
para permitir, ou não, tal deslizamento. Torna-se tam- 
bém razoável aceitar que os estados de tensão (interna) 
existentes na textura antes da carga e descarga influen- 
ciam esto comportamento. Em determinados casos e para 
certos ângulos de tensões máximas, o deslizamento não 
pode ocorrer e, neste caso, a deformabilidade é mínima. 

O modelo de deformação por deslizamento já foi 
sugerido para fissuras finas de Griffiths (º), ('º), mas pa- 
rece que se poderão vir a desenvolver teorias mais com- 
plexas. Espera-se que o conhecimento experimental aqui 
indicado contribua para estabelecer um modelo matemá- 
tico razoavelmente correcto para a deformabilidada de 
maciços rochosos 


pondentes de «E» num modelo isotrópico, isto é aos va- 
lores de «E» indicados no diagrama central da fig. 1 para 
za =0 e a = 90º, para os mesmos valores do parâmetro 
C (ver artigo seguinte). Os valores determinados de am- 
bas as formas são comparados no quadro seguinte. 

Estes resultados estão de acordo com a ideia que, 
para a fase de projecto (no caso de Las Portas), a reali- 
zação de ensaios com as medições necessárias para de- 
terminar os valores de E e E' (considerando a monotro- 
pia) irá, na prática, conduzir aos mesmos valores que os 
que implicam medições mais reduzidas e que podem ser 
interpretadas dentro da teoria mais simples da isotropia. 
Considera-se que esta conclusão pode ser extensível a 
todos os casos nos quais a anisotrraia é pouco signifi- 
cativa, seja qual for o modelo de anisotropia. 

Os ensaios, em laboratório, da resistência triaxial 
e da deformabilidade em arm.ostras de rocha matriz, ha- 
bitualmente considerados nos programas de estudo da 
deformabilidade dos maciços de fundação das barragens, 
têm, basicamente, o propósito de complementar os re- 
sultados dos ensaios «in situ». Foram apresentados al- 


Valores 
dos módulos de elasticidade (kg cm) 


Valores 
Formação n.º da E (mi) 
paso E 
[e] 19 76 000 
1 24 681 000 


E «E» (a = 0º) «E» (a =90º) 
55 000 70 000 52 000 
48 500 90 000 42 000 


2 — EFEITO DA DIRECÇÃO SOBRE A DEFORMABILI- 
DADE (OU «MÓDULO DE ELASTICIDADE») DE 
UM MACIÇO ROCHOSO E DO MATERIAL ROCHA 


A fig. 1 mostra resultados de ensaios em amostras 
em carga triaxial o de ensaios em galerias, em várias 
direcções, no local da barragem de Las Portas, em Es- 
panha, Deixando de lado os efeitos da pressão lateral 
em amostras em carga triaxial, é evidente que, sendo 
os «módulos de elasticidade», «E», em diversas direc- 
ções, de maciços rochosos, definidos por curvas do tipo 
indicado nessa figura, ou mesmo do tipo indicado na fig. 
4, para ensaios «in situ» ou em laboratório, a anisotropia 
está longe de ser do tipo mecnotrópico. Tal hipótese não 
vale, mesmo como aproximação teórica. Deve, além 
disso, notar-se que a anisotropia de «E» não é aquela que 
corresponde a curvas como as indicadas na fig. 1. De 
facto, os valores de «E» aí representados foram calcu- 
lados usando fórmulas válidas só para sólidos isotró- 
picos. 

No caso de Las Portas, foi demonstrado (*) que, 
para as duas formações 1 e 6, os valores E e E' do 
maciço rochoso, obtidos considerando um modelo mono- 
trópico (*), são praticamente iguais aos valores corres- 
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guns resultados nos Congressos da Associação Interna- 
cional das Rochas (*), (*). Uma das conclusões práticas 
mais importantes obtidas da comparação entre os resul- 
tados dos ensaios em laboratório e «in situy é que as 
caracteristicas de deformabilidade de um maciço ro- 
choso não podem ser avaliadas apenas através de en- 
saios de amostras em laboratório. Pelo contrário, podem 
ser cometidos erros importantes em relação aos valores 
dos módulos de elasticidada e da anisotropia da massa 
rochosa. À fig. 1 demonstra clãramente essa ideia. Esta 
figura compara, em representação polar, a deformabili- 
dade do maciço rochoso e da rocha matriz em amostras 
cilíndricas com 12 cm de altura e 6 cm de diâmetro en- 
saladas uniaxial e triaxialmente. À formação em questão 
é um horizonte de xistos no local da barragem de Las 
Portas. Duas conclusões podem ser tiradas imediata- 
mente dessa comparação: 
— os valores de «E» são menores para o maciço 
rochoso; 
— os dois tipos de ensaio mostram anisotropia di- 
ferente sendo mais acentuada nas amostras da 
rocha. 


Para amostras com planos de xistosidade paralelas 
ac eixo longitudinal (a = 0º), os valores de «E» são mais 
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baixos que os obtidos no caso de amostras com planos 
de xistosidade perpendicular aos seus eixos. Provavel- 
mente, isto é devido ao facto das amostras dos furos 
paralelos à xistosidade serem mais afectados pela ex- 
tracção das amostras. 


Um outro aspecto importante destes ensaios resulta 
du geometria das amostras e dos métodos de ensaio, 
que modificam a verdadeira anisotropia. De facto (ver 
figs. 1 e 2-1), considerando d o diâmetro da amostra 


e h a sua altura, o ângulo 5 é dado por 
d 
G=arc tg — 
h 


Da fig. 2 pode-se ver que para valores de q tais 
que O ac< à, não pode haver praticamente um contri- 
buto importante para a deformação longitudinal do desli- 
zamento ao longo dos planos de xistosidade passando 
através do centro do extensómetro, desde que as placas 
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De facto, nessas circunstâncias, o atrito por si só já não 
é suficiente para evitar o deslizamnto, Isto acontece no 
intervalo B das curvas da fig. 1, onde «E» é mínimo. Para 
valores de a, tais que 7/2 —- B < a< 90º, já não ha- 
verá teoricamente deslizamento ao longo dos planos 
de xistosidade, já que apenas o atrito será então sufi- 
ciente para o evitar. Na prática, o contributo para a de- 
formação do deslizamento tangencial ao longo dos planos 
da xistosidade diminuirá, dando azo a um aumento de 
«E», conforme se pode ver no intervalo C das curvas, 

No caso do maciço rochoso( diagrama central das 
figs. 1 e 2-Il) a situação é mais fácil de explicar. De 
facto, haverá sempre alguns planos de xistosidade aber- 
tos através dos quais pode haver deslizamento, uma vez 
quo as forças resistentes e o atrito forem vencidos. 
Assim, sempre que « for tal que O <a < m/2-8B, ha- 
verá deslizamento ao longo dos planos de xistosidade 
abertos. O número destes planos diminui conforme a se 
aproxima a 7/2 — &, de modo que o contributo desse 
deslizamento para o deslocamento medido no centro 
da superfícia em carga também diminui (intervalo D da 
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Fig |- Comparação entre anisotropia de deformabilidade 
de massas rochosas e do material rocha 


de carga nas bases das amostras evitem tal desliza- 
mento. É óbvio que este contributo, que é inexistente 
quando «a = O, aumentará na prática conforme a se apro- 
xima de 5, (intervalo À das cirvas da fig. 1 correspon- 
dente ao material rocha). Se o ângulo de atrito interno 
nos planos de xistosidade for À para valores de a, tais que 
ECalm/2-B | ecom7/2-A< 8), haverá teori- 
camente um contributo mais importante para a deforma- 
ção longitudinal através do deslizamento ao longo dos 
planos de xistosidade que atravessam os extensómetros. 


curva da fig. 1 referente ao maciço rochoso) e conse- 
quentemente o «módulo de elasticidade» aumenta. Para 
a=m./2- 8, teoricamente não haveria mais desliza- 
mento, sendo o mesmo verdade para valores de «a tais 
que 7/2-—- 8 <a < 90º. Por outro lado, contudo, o con- 
tributo, para o deslocamento, do fecho dos planos de es- 
tratificação e xistosidade abertos aumentará. O intervalo 
E da curva corresponde àquele intervalo de a, igual 
ou muito próximo de 7/2 — 5, o que de facto acontece 
no caso em estudo. 


[I-ENSAIOS 


|- ENSAIOS DE MATERIAL ROCHOSO 


INTERVALO (A) DE Q 
o<a<6 


IN SITU. 


Wi) 
Hi MI IH ; 


ROC ER 


INTERVALO DE Q 
o<a<T,p 


a =p 


b<a <p 


A 
: fil 


INTERVALO (0) DE O 
a B<a< 90" 


HU) | 


INTERVALO (E) DE Q 
Tp<a<90º 


Fig-2 Representação esquematica da influência do a na 


contribuição para a deformabilidade do escorregamento 


ao longo dos planos de xistosidade 


3 — INTERPRETAÇÃO DOS RESULTADOS DOS EN- 
SAIOS «IN SITU» CONSIDERANDO PARÂMETROS 
DE DIRECÇÃO E DE PROFUNDIDADE 


Para além de certas condições locais, tais como 
e proximidade de falhas ('), túneis (*), etc., a experiên- 
cia mostra que a profundidade, a orientação e a inclina- 
ção da superfície em carga e a posição da câmara de en- 


saio em relação à superfície do terreno sobrejacente 


têm uma influência fundamental nos valores de «E». 
Cada superfício carregada pode ser caracterizada por 
dois parâmetros geométricos. O primeiro é o ângulo « 
entre o normal à superfície em carga e os planos das 
principais discontinuidades do maciço rochoso (planos 
de xistosidade, de sedimentação ou de diaclases, em 
especial da família mais importante). O segundo é o 
parámetro C, representando a cobertura efectiva da cà- 
mara de ensaio. Este parâmetro é função das distâncias 
verticais e horizontais entre a câmara e a superfície do 
terreno. 
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Se for efectuado um número de ensaios suficientes 
em câmaras de orientação diversa, pode ser traçado um 
conjunto de curvas do «módulo de elasticidade» em fun- 
ção do parámetro C. Cada curva corresponde a um va- 
lor do parámetro a. Às curvas são obtidas a partir dos re- 
sultados dos ensaios e, em cada caso específico, são 
cada vez mais fáceis de traçar quando comparadas com 
as obtidas em estudos já realizados para outras forma- 
ções rochosas (foi só através de tais comparações que 
se tornou possível traçar as curvas para o local da bar- 
ragem de Janovas). À fig. 3 mostra as curvas obtidas 
para diversas formações rochosas ensaiadas em quatro 
locais de barragem, 

Para os casos estudados, o parâmetro geométrico 
ec foi determinado por várias fórmulas, tendo em conta 
a: posição da câmara de ensaio por baixo da superficie 
do terreno. Assim, nos locais de Las Portas e Sfikia foi 


2v + h 


achada a expressão e para o local de La 
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1,bv + h 
Baells a expressão ——— — . Em ambos os casos C é 


apenas uma média pesada entre as distâncias horizontal, 
h, e vertical, v, até à superfície do terreno, sendo a úl- 
tima considerada como um peso de 2 nos casos de 
Sfikia o Las Portas e de 1,E no caso de La Baells. Nos 
vales muito regulares, apenas com ligeira variação de 
secção, há uma certa relação entre os valores dos pesos 
que devem ser considerados para v e o gradiente da 
compressão natural em profundidade: quanto menor o 
último maior o primeiro. Nesse aspecto, Janovas é um 
caso algo especial. A cobertura h é de facto considerada 
como tendo um peso que é uma função de h' (cobertura 
horizontal perpendicular a h), com um limite superior de 
1, Quer isso dizer que o peso de h foi sempre, neste 
caso, menor que ou igual a 1. Este peso tenta considerar 
a influência da discontinuidade frontal do horizonte cal- 
cário 1 sobre a anisotropia dos diversos horizontes. As 
curvas relativas a Janovas, que se vêm na fig. 3, indicam 
muito claramente essa influência. De facto, enquanto os 
horizontes 3 e 5, longe de 1, mostram uma maior defor- 
mabilidade quando a carga é paralela à estratificação, 
o horizonte 2 tem comportamento quase isotrópico. É 


evidente, também, uma transição progressiva da aniso- 


tropia do horizonte 1 para a dos horizontes 3 e 5. 

Os conjuntos de curvas do tipo das que aparecem 
na fig. 3 são de bastante interesse prático, uma vez que 
permitem traçar os diagramas de «E», «E», etc., ao 
longo da fundação, para várias direcções ou apenas na 
direcção da normal à superficie de fundação projectada, 
«E, ». Para este efeito, basta escolher um número apro- 
priado de pontos para essa superfícia e calcular para 
cada um: 


1— o valor « que corresponde à normal à super- 
fícia da fundação num dado ponto; 


2 — os valores de v o h, desprezando a presença 
de qualquer recobrimento anómalo local, ou 
afloramento; 


3 — o valor do parâmetro C e; 


4 — o valor de «E» que corresponde aos valores de 
C e « obtidos das curvas da fig. 3. 


A fig. 4 mostra as curvas dos módulos de elastici- 
dade E, » ão longo da superfície de fundação de di- 
versos locais de barragens em Espanha e na Grécia. 


4 — IMPORTÂNCIA DA INJECÇÃO NA DEFORMABI- 
LIDADE 


Considera-se geralmente que a injecção de caldas 
de cimento melhora as propriedades dos maciços rocho- 
sos ('), mas poucos resultados têm sido apresentados. 
A fig. 5 mostra uma comparação entre a anisotropia da 
deformabilidade antes e depois do tratamento por injec- 
ção, efectuado no local da barragem da Sfikia, em an- 
fibolitos xistosos. Este resultado pode talvez ser consi- 
derado típico de maciços rochosos em camadas. De 
facto, na barragem de Cambambe (Angola), foram obti- 
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dos resultados idênticos, não publicados, em arenitos em 
em camadas. De uma análise dos resultados apresenta- 
dos na fig. 5, pode-se deduzir o seguinte: 


— Os valores de «E» relativos à superficie normal 
à xistosidade são substancialmente mais altos 
do que os obtidos para superfícies paralelas à 
xistosidade (carga perpendicular aos planos de 
xistosidade). 


— O valor minimo de «E» antes da injecção cor- 
responde a um valor de « da ordem dos 50º. 
Esta valor do ângulo a representa a situação 
para a qual o contributo do deslizamento ao 
longo dos planos de xistosidade (ou microfis- 
suras) é máximo. Deve-se frisar que, apenas 
quando 7/2 — a é maior que A (sendo 4 o ân- 
gulo de atrito interno nos planos de xistosidade), 
haverá deslocamento por deslizamento. Isto in- 
dica que o ângulo A, sob condições normais de 
humidade, será da ordem dos 40”. De facto, 
em ensaios de corte no local, foi obtido um 
valor de 32º com o maciço rochoso totalmente 
saturado. 


— O tratamento de injecção contribui efectivamente 
para melhorar o «módulo de elasticidade» do ma- 
ciço rochoso. Essa melhoria é reflectida nos 
ensaios no local de Sfikia através de um au- 
mento do «E» de cerca de 15 a 25 % para su- 
perfícies paralelas à xistosidade (a — 90º); cerca 
de 100 % para superfícies normais à xistosidade 
a=>0)e110sa115% para superfícies corres- 
pondentes a «= 33 a a = 55º, respectivamente, 


— O tratamento do local por injecção acentuou a 
anisotropia. 
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“a  MURVA DE É tantes da injecção! 


5 — DEFORMAÇÕES IRRECUPERÁVEIS DOS MACIÇOS 
ROCHOSOS 


Conforme mencionado anteriormente, durante os 
ciclos de carga e descarga dos ensaios de deformabili- 
dade observa-se sempre uma certa deformação irrecupe- 
rável, Este facto tem sido atribuido ao fecho das dia- 
clases e microfissuras que se mantêm parcialmente fe- 
chadas durante a descarga. Tem sido indicado (") que 
esta deformação ou deslocamento permanente (assen- 
tamento) é reduzido pela injecção e, por outro lado, sa- 
be-se que, na prática corrente, a injecção enche apenas 
as diaclases acima duma certa abertura (0,1 a 0,2 mm). 

Os deslocamentos permanentes no fim dos ciclos 
de carga/descarga dão, normalmente, uma indicação da 
de compressão dos maciços rochosos. À ideia de estudar 
c anisotropia dos deslocamentos totais permanentes teve 
concretização pela primeira vez, no caso do local de 
Sfikia. Tentou-se uma correlação entre a média dos des- 
locamentos totais permanentes de cada par de superfl- 
cies paralelas em carga e o parámetro C (fig. 6-1). 
Essa figura mostra que: Eae 


DESLOCAMENTO imm) 
ho 


10 15 20 25 30 
PARAMETRO Caim 


NOTA 


-Os valores de M e C 
constam da Fig.3 


mostra que a anisotropia dos deslocamentos permanentes 
está relacionada com a deformabilidade: 


— Os mais altos valores de «E» correspondem aos 
deslocamentos permanentes mais pequenos e 
vice-versa, 


— Os deslocamentos permanentes maiores ocorrem 
em relação a uma gama de valores a para a qual 
o potencial de deslizamento ao longo dos planos 
de xistosidade é máximo. 


— À injecção reduz os deslocamentos permanentes. 


6 — COMPORTAMENTO OBSERVADO EM FUNDAÇÕES 
DE BARRAGENS 


O comportamento real dos maciços rochosos das 
fundações das barragens, quer durante a construção (de- 
vido à sobrecarga com o peso da obra), quer em serviço 
(devido às cargas transmitidas pela barragem e também 
às pressões de água intersticial), é actualmente obser- 


Curva de deslocamento 
permanente total 


Fig.6-Resultados obtidos no local da barragem de Sfikia. 
(I) Correlação entre o deslocamento permanente total 


e o parametro CIT) Anisotropia do deslocamento 
permanente em planos normais à xistosidade. 


— Os deslocamentos permanentes diminuem com 
o aumento de C€C. Os dados indicam que o deslo- 
camento total permanente é reduzido em mais 
que metade quando C aumenta de 15 para 25 m. 


— O deslocamento permanente apresenta um com- 
portamento anisotrópico que é bem definido pelo 
parâmetro a. 

O diagrama polar da fig. 6-|l foi obtido para um 

valor C igual a 20 m. Uma comparação com o da fig. 5 
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vado com diversos tipos de instrumentos ("2), 

Além das medidas de deslocamentos (usando pen- 
dulos invertidos), inclinações (clinómetros) e desliza- 
mentos (medidores de deslizamentos), foi dada especial 
atenção às medições de deformações (extensões). Essas 
medições são feitas com extensómetros eléctricos com 2 
ou 3m de comprimento, idênticos aos extensómetros 
normais para betão ("'). Os autores já descreveram ("), 
(*) essa aparelhagem e as normas de formação de gru- 
pos bi-dimensionais, de quatro extensómetros; de grupos 
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tridimensionais: da extensúmetros compensadores (para 
efectuar medições das mudanças de comprimento dos 
maciços rochosos não relacionadas com tensões, mas de- 
vidas à temperatura e a mudanças de comprimento «hi- 
grométricas», bem como as técnicas apropriadas para 
colocar esses extensómetros nas fundações. 

Algumas das conclusões mais interessantes anterior- 
mente obtidas das medições de extensões em maciços 
rochosos foram (") (1): 


— q acordo, em ordam de grandeza, dos «módulos 
de elasticidade» obtidos pelos ensaios «in situ» 
e dos calculados a partir das extensões medidas 
nas fundações rochosas; 


— os aumentos do «módulo de elasticidade» verti- 
cal com a carga, observados durante a constru- 
ção de algumas barragens que foram providas 
de aparelhagem; 


-—— as grandes deformações dos maciços rochosos 
durante os trabalhos de injecção, conforme ob- 
servado em várias barragens. 

Agora podem ser acrescentados mais alguns resul- 

tados e conclusões, tais como: 


— dados interessantes sobre a anisotropia da defor- 
mabilidade das fundações de algumas barragens; 


— em geral, uma grande influência das mudanças 
de comprimento «higrométricas» do maciço ro- 
choso (devido a variações de humidade nos po- 
ros ou à penetração de água na textura da ro- 
cha) sobre os seus estados de tensão; 


— a existência da um atrazo no tempo entre as si- 
tuações limita das acções mecânicas (peso má- 
ximo ou fim de construção, níveis mínimos ou 
máximos da albufeira) e o limite correspondente 
dos valores (máximo ou minimo) das extensões 
no maciço rochoso; 


— uma redução no «módulo de elasticidade» com a 
carga para as direcções da tensão de corte má- 
xima nas quais o deslizamento (ou as microfis- 
suras) nos planos de xistosidadea contribui clara- 
mente para a deformabilidade. Este aspecto 
parece ser influanciado pela maior ou menor 
restrição imposta a esse deslizamento pela geo- 
metria dos blocos da barragem. 

Nos pontos onde os estados de tensão nas funda- 
ções são facilmente calculados (por exemplo as bases 
dos blocos de betão) e onde existem grupos de extensó- 
metros, foram obtidos os valores dos «módulos de elas- 
ticidade» em várias direcções, o que permitiu uma visão 
do carácter anisotrópico da rocha de fundação. Isto foi 
tornado possível por uma comparação dos estados de 
tensão calculados e os estados de tensão indicados pelos 
extensómetros. No caso de barragens abóbada, os pontos 
da fundação normalmenta usados para este efeito foram 
os pontos do bloco central no plano vertical axial, en- 
quanto que para barragens de gravidade estes pontos 
foram geralmente os do plano vertical axial dos blocos 
mais altos. Observou-se que, em diversas instâncias, os 
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valores dos «módulos de elasticidade» obtidos ao longo 
das diferentes direcções, correspondem perfeitamente 
à maior ou menor deformabilidade prevista de acordo 
com as atitudes dos planos de discontinuidade (xistosi- 
dade, estratificação, diaclases, microfissuras ou seme- 
lhantes). Um bom exemplo é o dos «módulos de elasti- 
cidade» obtidos na fundação da barragem de El Atazar 
(Espanha). Para esse fim, o estado de tensão criado 
pelo peso do bloco central foi comparado com o estado 
de tensão correspondente medido por extensómetros. 
Foi calculado um «estado de tensão», considerando um 
sólido continuo, homogéneo e isotrópico, carregado su- 
perficialmente ao longo de uma faixa com comprimento 
indefinido e uma largura constante e igual à largura mé- 
dia do bloco central, A rocha da fundação da barragem 
de El Atazar é um xisto com planos de xistosidade quase 
verticais e uma direcção praticamente perpendicular ao 
eixo do rio. Foram obtidos valores de 60 000 kg. cm —< 
E de 40 000 kg. cem—2 para os «módulos de elasticidade» 
em direcções paralela e normal à xistosidade, respecti- 
vamente. 

As mudanças de comprimento «higrométricas» são 
facilmente detectadas, como se disse, por extensómetros 
compensadores colocados dentro dos maciços rochosos 
das fundações. Elas dependem, obviamente, da natureza 
da rocha. Tais mudanças de comprimento são geral- 
mente pequenas em rochas sãs, aumentando com o grau 
de alteração, Para rochas ligeiramente alteradas, a ex- 
pansão «higrométricar devida à embebição e saturação 
da rocha (ou, por outras palavras, a uma expansão au- 
togénia devida à presença de água) podes atingir valores 
d3 + 200.106 a + 300.10 —€ [u/m). Esta expansão 
dc maciço rochoso parece ser benéfica nara o compor- 
tamento da barragem. De facto, contribui para a redu- 
ção da permeabilidade do maciço rochoso de fundação, 
diminuindo a filtração através dele, produz auto-com- 
pressão, dando azo a uma redução e mesmo à eliminação 
de possíveis tensões de tracção no pé (base de mon- 
tante) não só no maciço rochoso, mas também na parte 
inferior dos blocos da barragem e cria uma certa com- 
pressão no betão perto das fundações, Ou seja, no caso 
das barragens abóbada, no soco da base da abóbada. 

O intervalo de tempo, ou atrazo, entre as situações 
limite das acções mecânicas e as de deformações no ma- 
ciço rochoso tarda, por vezes, algumas semanas. Pensa- 
-se que está relacionado com a natureza plástica e visco- 
elástica das deformações sofridas pelo maciço rochoso, 
especialmente as que resultam das discontinuidades 
imicrofissuras, diaclases finas, caixas de falhas, etc.). 
É evidente que quanto mais denso o programa de leituras 
durante estes períodos, maior é a precisão com que este 
atraso pode ser determinado, 


A fig. ? mostra alguns resultados das extensões 
obtidas por rochas das fundações de três barragens 
abóbada durante o primeiro enchimento das suas albu- 
feiras. Referem-se (') a uma pequena barragem-abóbada 
[EI Vellon, Espanha), (*) a uma barragem abóbada de ta- 
manho médio (Funil, Brasil) e (') a uma outra muito 
alta 6 com um coroamento muito comprido (Almendra, 
Espanha). As rochas da fundação são respectivamente, 
gneiss micácio, gneiss e granito. | 
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Na fig. 6 da referência ("), já tinham sido indica- 
dos, para o caso da barragem de El Vellón, resultados 
idênticos, embora referentes a uma fase anterior do en- 
chimento da albufeira. As observações feitas naquele tra- 
balho não são, por isso, repetidas aqui, sendo agora indi- 
cadas apenas certas conclusões complementares. 

No grupo situado perto da base de montante 
do primeiro bloco da margem direita, notou-se que a de- 
formação na direcção vertical aumentava consideravel- 
mente em relação à situação indicada na fig. 6 de ("). 
Isto foi devido ao efeito da injecção de cimento na fun- 
dação daquela zona da parte inferior da encosta. No 
caso dos dois grupos simétricos situados a nível médio 
das encostas perto do pé de jusante da abóbada, notou-se 
que, em geral, as deformações aumentavam em propor- 
ção com o aumento da carga hidrostática. As maiores 
deformações de compressão continuam a ocorrer ao 
longo da direcção no plano vertical radial que baixa a 
45º para jusante (direcção 5) e das direcções normais 
às vertentes no plano tangencial (direcção 6 na margem 
esquerda e direcção 7 na margem direita). 

No caso da barragem do Funil, deve-se notar que, 
como os grupos da fundação desta barragem não têm 
extensômetros compensadores, não é possível corrigir 
os valores dessas deformações pela redução das mudan- 
cas do comprimento de origem termo-higrométrica, Isto 
não permite uma avaliação das deformações devidas a 
tensões (extensões). As mudanças de deformação indi- 
cadas na fig. 7 para esta barragem são, de facto, mudan- 
cas de comprimento total unitário, isto é, mudanças de 
comprimento medidas, divididas pelo comprimento do 
extensómetro (extensões, mais expansões «termo-higro- 
métricas»). A partir do seu estudo cuidadoso, pode-se 
concluir o seguinte: 


1 —— Todas as mudanças de comprimento totais unitárias 
são positivas. Por conseguinte mostram claramente 
uma expansão da rocha por saturação. 


2 — No único grupo existente a jusante (pé do bloco 
central), as mudanças de comprimento unitárias 
positivas indicadas na figura devem ser mais baixas 
do que a expansão efectiva, especialmente para 
aquelas direcções mais influenciadas pela compres- 
são. Pelo contrário, em relação aos grupos de mon- 
tante, especialmente para as direcções onde ocorre 
maior descompressão, essa descompressão deve 
aumentar as mudanças de comprimento unitárias 
positivas. De facto, o valor médio das mudanças 
de comprimento unitárias reais indicadas pelos 
extensômetros do grupo de jusante é de + 4.10-6, 
enquanto o valor médio das mudanças de compri- 
mento unitárias dos extensómetros dos grupos a 
montante varia de um minimo de + 45.10 —6 a 
um máximo de + 120.10-6, 


3 — Nota-se claramente que a expansão «higrométrica» 
é mais acentuada na fundação da margem esquerda 
do que na da margem direita. Isto tem uma relação 
estreita com a mais fácil percolação de água na 
margem esquerda, causada pelo maior número e 
maior importância dos acidentes tectónicos que a 


afectam. Considerando os resultados das medições, 

parece também que é mais elevada a expansão 

«higrométricar a cotas mais altas em ambas as 

margens. 

Para a barragem de Almendra, uma análise dos 
dados indicados na fig. 7 permite tirar as seguintes con- 
Clusões preliminares: 


1 — No bloco central, a montante, as extensões positivas 
indicam a existência de tensões de tracção na fun- 
dação devidas à carga máxima da água na barra- 
gem. Estas extensões foram o resultado de uma 
descompressão que tinha ocorrido na rocha ao longo 
de um periodo de um ano a partir de 24/4/72 
(albufeira à cota 720,4 m). As extensões positivas 
principais ocorrem numa direcção quase vertical, 
com inclinação para montante. Dentro deste periodo 
dao um ano, houve dois sub-períodos, de dois meses 
cada, durante os quais começaram a maior parte 
das deformações. Estes sub-períiodos coincidem com 
e princípio de uma subida mais rápida na albufeira. 
Houve também aumantos repentinos nas pressões 
de águas intersticial na mesma altura na área das 
fundações dos blocos centrais. 


? — Nos blocos centrais, a jusante, foi precisamente du- 
ranto os subperiodas indicados que as extensões de 
compressão se desenvolveram mais rapidamente, 
O aumento de carga transmitido pela abóbada 
aquela área da fundação tornou-se assim óbvio. As 
deformações de compressão são mais altas ao longo 
das direcções 3 e 5, inclinando para jusante. 


*-— Nos grupos tridimensionais situados na parte de 
jusante da fundação de alguns blocos laterais e para 
além das direcções 3 (vertical) e 1 (tangencial, ho- 
rizontal), são também privilegiadas as seguintes 
direcções em relação às extensões desenvolvidas: 
na margem esquerda, direcção 6 (aproximadamente 
normal à superfície da fundação) e direcção 9 (ho- 
rizontal a jusante); na margem direita 7 e 8, quer 
dizer, as simétricas das direcções 6 e 9 respectiva- 
mente. ds 


£ — Na direcção radial 2, as extensões são, em geral, de 
menor valor. 
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Um modelo gaussiano para a dispersão 
de poluentes na atmosfera 


RESUMO 


Apresenta-se um modelo Gaussiano de dispersão 
de poluentes na atmosfera, considerando o efeito de 
fontes múltiplas pontuais, a existência de inversão de 
temperatura em altitude, a variabilidade da direcção e 
velocidade do vento e da estabilidade da atmosfera, 

4 estabilidade da atmosfera é classificada segundo 
& método de Pasquill sendo o desvio padrão da concen- 
tração de poluentes calculado a partir de aproximações 
numéricas de gráficos de Pasquill-Giftord. 

No final são apresentados alguns exemplos que 
ilustram o efeito da variação da estabilidade da atmos- 
fera, da existência de inversão e a interferência de 
fontes próximas. 


O. INTRODUÇÃO 


0.1. Ao considerar aspectos estritamente técnicos da 
poluição do ar ou da água, duas vias complemen- 
tares imediatamente se nos oferecem: a da geração 
dos poluentes e a da sua dispersão. 

Nesta comunicação restringir-nos-emos a uma 
faceta muito particular da dispersão de poluentes 
na atmosfera que, por consistir num método de 
cálculo expedito, se presume possa corresponder 
aos interesses práticos dos participantes neste con- 
gresso. 

Métodos deste tipo, embora mais simplifica- 
dos, são hoje dos mais correntes e neles se baseiam 
praticamente todos os códigos de dimensiona- 
mento de chaminés actualmente existentes. Tais 
códigos consideram em geral uma fonte única e 
uma direcção de vento constante. O modelo aqui 
adoptado considera o efeito de fontes múltiplas, 
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ABSTRACT 


This report presents a Gaussian dispersion mode! 
for gaseous pollutants which accounts for the effects 
of multiple ponctual sources, temperature inversion, at- 
mospheric stability and variations in both wind speed 
and direction. 

The stability of the atmosphere is classified into 
Pasquill's classes and the standard deviations ot the 
concentration ( 90 jare calculated bv means of a 
numerical approximation to Pasquill-Gifford's plots. 

A few examples ilustrating the effect of atmos- 
pheric stability and of inversion on ground level concen- 
tration and of the interference of nearby sources are also 
worked out. 


a variabilidade da direcção do vento, as classes 
de estabilidade da atmosfera e o efeito da altura 
de inversão. 

Relativamente à validade da sua aplicação e à 
justificação do seu uso remete-se para [1] onde se 
discutem as bases de códigos actualmente exis- 
tantes. 


1. MODELOS GAUSSIANOS 


1.0. Do modo a situar os modelos Gaussianos na sua 
perspectiva, resumem-se seguidamente as considera- 
ções fundamentais que levaram à sua formulação. 
Basicamente, tais considerações são: complexidade 
excessiva das equações fundamentais que regem 
os fenómenos de difusão; ausência de conheci- 
mento suficiente quanto à difusão turbulenta; im- 
precisão quanto a dados meteorológicos e caracte- 


303 


Tels 


rísticas das fontes de emissão; custo excessivo das 
simulações com modelos mais complexos nara sua 
utilização sistemática. 

Todos os aspectos referidos têm sido e são ob- 
jecto de trabalho no grupo de investigação onde 
nos integramos. A complexidade das equações não 
constitui limitação essencial face às técnicas numeé- 
ricas entretanto desenvolvidas, apenas o consti- 
tuindo os computadores de que podemos dispor. A 
investigação sobre a difusão turbulenta, embora não 
tenha ainda permitido um completo esclarecimento 
do assunto, permite, no entanto, em conjugação com 
as técnicas numéricas desenvolvidas, a solução ade- 
quada duma classe cada vez maior de problemas. 
As limitações quanto a dados meteorológicos, so- 
bretudo de natureza local irão sempre existir; toda- 
via, a conjugação dos conhecimentos relativos à 
turbulência com as técnicas numéricas permitiram 
já um avanço considsrável na previsão de micro- 
climas e da sua interacção com as alterações locais 
motivadas pelo relevo, a florestação, os poluentes 
introduzidos, etc. Permanecem, pois, como muito 
importantes, as limitações relativas ao custo das 
simulações em computador, as que se referem à 
inventariação das fontes e a informação meteoro- 
lógica relativa às médias e grandes escalas, 

Todos os aspectos referidos são hoje objecto 
de investigação em centros importantes. Entre nós, 
a eles se ligam as publicações [2] a [9]. 


Com os modelos Gaussianos procura-se obter a 
universalidade e segurança que conferem as equa- 
ções fundamentais, evitando simultaneamente a 
complexidade e custo que resultaria da sua solução 
tendo em conta todos os factores, Com essa fina- 
lidade, as equações fundamentais são simplificadas 
admitindo: 


— regime estacionário 
— coeficientes da difusão constantes 
— vento uniforme e de direcção constante 


— ausência de reacção quimica entre po- 
luentes 


— difusão desprezável na direcção do vento 


Com estas hipóteses, a equação fundamental 
é linear e torna-se possivel obter soluções analíti- 
cas simples. Dadas as simplificações introduzidas, 
tais soluções só muito grosseiramente descrevem 
o fenómeno, Para que os resultados se aproximem 
da realidade, a via seguida é a de ajustar empiri- 
camente a solução analítica aos resultados experi- 
mentais através da dependência dos coeficientes 
de difusão das coordenadas espaciais. De acentuar 
que tal dependência é introduzida na solução e não 
na equação de partida. Às razões deste modo de 
proceder são essencialmente de natureza prática, 
resultando a sua utilidade da validação experimen- 
tal que foi possível conferir-lhe. 

É sobretudo de realçar que, embora de forma 
aproximada, é possivel ter em conta nesta formula- 


ção a influência da estabilidade da atmosfera, a 
qual se traduz na forma da dependência espacial 
dos coeficientes de difusão, Tal denendência foi 
inicialmente introduzida por Pasquill e posterior- 
mente modificada por Gifford. A simplicidade do 
tratamento e a extensa validação experimental a 
que foi submetida sobretudo nos E.U.A. e no 
Reino Unido estão na origem da extensa utilização 
que delas hoje se faz 


1.2. Apesar da validação experimental e da sua gene- 
ralidade, os modelos Gaussianos apresentam intrin- 
secas limitações, que são frequentemente esqueci- 
das. A utilização duma simulação numérica mais 
elaborada, tal como a própria verificação experi- 
mental revelam facilmente que o relevo, o perfil 
do vento, as características do terreno (le da vege- 
tação) etc. introduzem desvios que o modelo não 
pode explicar, 

Tal não impede, porém, a sua utilidade em es- 
timativas preliminares pelo que se justifica não só 
o seu uso generalizado como a sua extensão a si- 
tuações não previamente consideradas. 


2 MODELO ADOPTADO 


2.1. Tendo em conta as limitações mas reconhecendo 
as vantagens da formulação de tipo Gaussiano, 
o caso de fontes múltiplas e da variabilidade da 
direcção do vento foi considerado. Além disso, é 
feita a sua extensão ao caso em que existe uma in- 
versão de temperatura. 

Tendo em conta as hipóteses típicas do mo- 
delo Gaussiano, pode obter-se para uma distribui- 
ção continua de fontes pontuais a seguinte solução 
analítica da equação fundamental para o caso em 
que no solo não existe absorção de poluente: 


q mn OQ, 
tt “É fe 4a) = — E E Rs Ae gia 
= Cybexç) Calxx) 
tz — py tz +h,) ly— y;y 
Xp — ————— a — 4-BXP — n * ex —— 
| ' o á 2 0; : 20º 


em que ox coincide com a direcção média do vento; 
oz com a direcção vertical e oy com a direcção 
transversal, sendo: 


U —— velocidade média do vento 


— potência da fonte | (mg/s) localizada no 


x,» Y; —— coordenadas da fonte i no plano oxy 
5; 9%, — desvio standard da distribuição do poluente 
Ea na direcção horizontal e vertical respecti- 
vamente 
h. —— altura efectiva da tonte i 


À coordenada h; coincide com a altura da 
fonte no caso de não haver elevação sensivel do 
poluente devido à sua velocidade de emissão ou à 
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sua temperatura. No caso geral de uma chaminé 
existe uma elevação Ah, pelo que, nesse caso, 
h, =h, + 4h, sendo h, a altura geomé- 
trica de emissão: A introdinddn é 'de h;, eo seu valor 


é uma aproximação que se discute aim [1]. 


2.2. No caso de existir uma inversão de temperatura à 
altura H > h, . o seu efeito pode ser aproximado 
por um impedimento à passagem dos poluentes para 
além de H, A solução analítica para este caso ob- 
tém-se pela técnica das imagens, considerando re- 
flexões sucessivas no plano do solo e no plano 
à altura H. Embora teoricamente se deva conside- 
rar um número infinito de fontes virtuais, verifica- 
-sa que a precisão é suficiente considerando apenas 
três pelo que a solução final é, para a concentra- 
ção ao nível do solo: 


o mn 20; iai 
C (xy 0)=—=2 ed 
» T Uj= 1 “y 0 á 9? 


hi ah y 
o [sã Jem [Si]. 
Co 4(2H-h927 | 
cof 527] 


2 Cr 
2.3. As soluções analíticas anteriores são válidas para 
o caso dos coeficientes de difusão serem uniformes 
como já referido; o que caracteriza os modelos 
Gaussianos é a introdução duma dependência fun- 


cional entre os desvios standard (ou «mutatis mu- 


tandis» dos coeficientes de difusão) e as coorde- 
nadas espaciais. Tal dependência, introduzida em- 
piricamente por Pasquill e modificada por Gifford 
tem em conta a turbulência atmosférica a qual é 
indirectamente introduzida através das classes de 
estabilidade atmosférica. Tal dependência é ex- 
pressa por aqueles autores sob forma gráfica das 
quais obtivemos uma aproximação analítica através 
de um ajustamento por mínimos quadrados de ex- 
pressões do tipo: 


log q = E a, x) 


As mesmas classes de estabilidade estão li- 


gadas com o desvio standard da direcção do vento 
para tempos de amostragem curtos (de 10 a 60 
minutos) e a seguinte relação empírica foi estabe- 
lecida por Gifford [23]: 


QUADRO II 
A Classe de 
Estabilidade 
(Graus) de Pasquill 
25 A 
20 5 
15 Ç 
10 D 
5 E 
2,5 F 


Nos quadros seguintes apresentam-se as clas- 
ses de estabilidade em função de parâmetros me- 
teorológicos e para cada classe os coeficientes a 

correspondente a q 8 9, 


QUADRO Ill 
log à SBg ta, log x+a, (log x) 
Classe de É a 
estabili. ão + 
A - .B8685 1,21747 
B - 94252 | 1.15522 
E - 93453 1.06553 
D o - 1.15973 1.09563 
go - 1.30518 1.11223 
' F - 1.55707 1.15887 


QUADRO | 


Velocidade Insolação durante o dia 
do vento 
m/s Forte Moderada 
< l A A — B 
2 A —B B 
4 E— E 
6 C C— D 
> 6 Cc D 


Condições nocturnas 
cobertura de nuvens 


Ligeira 
B > 4/8 « 3/8 
C E F 
C D E 
D D D 
D D D 
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2.4. 


3.0. 


e 


3.2. 


QUADRO IV 


o” aj 


log x + a: (log x)" + a; (log x)' 


Classe de 
estabilid. 


5.52 


- 3.59404 


B - 2.88414 3.66 
Cc same [as 
D - 1.97892 1.74 

E né - 1.88705 1.47 
206895 | 1551 


As expressões anteriores são aplicáveis quando o 
vento mantém uma direcção uniforme. Se a varia- 
ção de direcção é lenta, os efeitos cumulativos 
resultantes são pequenos e as expressões são uti- 
lizáveis a cada instante desde que se introduza um 
sistema de coordenadas em que ox coincida sempre 
com a direcção do vento. Em relação a este sis- 
tema de coordenadas, a posição das fontes vai na- 
turalmente variar. Esta transformação é feita auto- 
maticamente na implementação do modelo em 
computador. 


IMPLEMENTAÇÃO DO MODELO EM 
COMPUTADOR 


O modelo foi implementado num calculador Hew- 
lett-Packard Modelo 30 com 5856 palavras de me- 
mória e equipado de um «plotter». A linguagem 
utilizada é o BASIC. 

Para aquele calculador, o número máximo de 
pontos da malha é 1296 e o de fontes não tem li- 
mite, provocando o seu aumento um acréscimo de 
tempo de calculador. A implementação do modelo 
noutra linguagem é, aliás, muito fácil. 

O programa calcula as concentrações em todos 
os pontos da malha, 

O traçado das isopletas é feito usando uma 
subrotina concebida para esse fim [211]. 

Os dados necessários para utilização do programa 
são: 


a) Malha: número de pontos segundo ox e oy e seu 
intervalo; 


b) Dados meteorológicos: classe de estabilidade, 
velocidade e rumo od vento e altura de inver- 
são, em função do tempo; 


c) Fontes: número de fontes; potência efectiva, ve- 
locidade de emissão, temperatura, localização. 


815 


915 


929 
620 


860 


700 


- 2.43013 


- 1.22382 


- 07458 


- 23946 01499 
- 13449 00076 
- . 15818 00306 


3.3. Na Fig. 1 apresentam-se os resultados parciais 


4.0. 


duma simulação na qual se utilizaram os dados 
seguintes: 

Fonte de altura efectiva 979m emitindo 
350 000 mgSO;/s. Não existe camada de inversão. 


Curvas 1 — vento NW; 2 m/s; classe de esta- 
bilidade A (muito instável). 


Curvas 2 — vento W; 3.5 m/s; classe de es- 
tabilidade B (instável). 

Traçaram-se as linhas de isoconcentração de 
SO. ao nível do solo. 


A concentração no solo é superior para a 
atmosfera mais instável, 

Fig. 2 —— A mesma fonte com vento de W 3.5 
m/s e classe de estabilidade B. Existe uma camada 
de inversão em altitude cuja base se situa a 600 m. 
A altura efectiva da fonte neste caso baixa para os 
580 m pois não pode ultrapassar a camada de inver- 
são. A concentração no solo é superior quando 
existo inversão em altitude. 


Fig. 3 — Consideram-se duas fontes de carac- 
terísticas iguais à da Fig. 1; vento W, 3.5 m/s; 


classe B. Não existe inversão. As fontes inter- 
actuam. 
CONCLUSÕES 


O modelo construido é simples de utilizar e 
versátil, apresentando em relação aos modelos 
Gaussianos comumente utilizados algumas exten- 
sões do maior interesse, nomeadamente a de ter 
em conta a altura de inversão. 

Este modelo, se pode considerar-se, no mo- 
mento actual, como um aceitável compromisso entre 
rigor, flexibilidade e custo, apresenta, como todos 
os modelos Gaussianos. limitações que importa não 
esquecer na consideração de situações concretas. 
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No modelo apresentado apenas se consideram 
fontes pontuais. A sua extensão a fontes distribui- 
das segue a mesma metodologia e não oferece 
dificuldades adicionais. Em todos os casos em que 
o rigor fornecido por este modelo é insuficiente, 
hã que recorrer aos métodos mais complexos e 
onerosos a que se fez referência. 
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FIG. 1 


CLASSE A VENTO W 2M/S 
CURVA 1 —.1 MG SO,/M: 


2 — 05 
3 — 01 
4 — .001 


vo pah 


CURVAS 
CURVAS 
& — 04 
BE — 02 


[E 


386 


FIG. 2 


LLASSE B VENTO W 3.5 M/5 
CURVAS 1 — INVERSÃO À 600M 
CURVAS 2 —S/ INVERSÃO 

A — 04 ME 50,/M! 

EP — .01 
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1— CLASSE A VENTO NW | 2M/S 
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”Cepstrum” (') 


Estudo experimental da sua aplicação ao 
sinal falado. ” Vocoders ”. 


RESUMO 


Este trabalho, que constitui relatório de um estágio 
realizado no Laboratório de Análise de Sinais do 1. A. C., 
é iniciado com uma análise resumida do sinal falado 
e um breve estudo teórico da técnica do cepstrum. Segui- 
damente é apresentada a aplicação desta técnica na de- 
terminação de alguns parâmetros de sinais daquele tipo. 
Por último, faz-se uma introdução ao estudo dos vo- 
coders. 


1 — INTRODUÇÃO 


Tendo conhecimento de um trabalho dos colegas 
Borges de Almeida e Silva Arroz que constitui o seu re- 
latório de estágio, trabalho esse que seguimos com inte- 
resse nalgumas das suas fases, sobre as técnicas de 
«Cepstrum», resolvemos dar continuidade ao estudo deste 
assunto, tanto mais que a justificá-lo estava o facto de 
ele ser recente e pouco divulgado, 

No trabalho referido, dos colegas Borges de Al- 
meida e Silva Arroz, faz-se um estudo da técnica de 
«Cepstrum» especialmente na sua aplicação à detecção 
de ecos e dos seus atrasos relativamente ao sinal 
«principal» e vários outros problemas com estes rela- 
cionados. Baseados nesse trabalho, na bibliografia nele 
indicada e ainda noutra que nos foi sugerida pelo Dr. Re- 
belo Simões resolvemos estudar o «Cepstrum» e as suas 
aplicações no domínio da fala, aquele em que a técnica 
de «Cepstrum» se mostrou mais frutuosa e onde teríamos 
algumas possibilidades dado o material que era posto 


(1) Leia-se «Kepstrum» 
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SYNOPSIS 


This paper, reprinted from a technical report of La- 
boratório de Análise de Sinais, 1. A.C., starts with a 
short description of speech signals and a brief theoretical 
analysis of cepstrum. These techniques are then applied 
to the determination of speech parameters. The paper 
ends with an introduction to vocoders. 


à nossa disposição pelo Laboratório de Análise de Sinais 
e que noutra parte indicamos resumidamente, 

Para tal estudo resolvemos adoptar o seguinte es- 
quema: 

No segundo parágrafo começamos por analizar o 
sinal falado, mostrando alguns exemplos e em seguida, 
no terceiro parágrafo, põe-se em evidência a necessidade 
o utilidade do «Cepstrum» no seu estudo. 

No quarto parágrafo estuda-se o «Cepstrum» em por- 
menor e as suas várias formas possíveis são apresenta- 
das e discutidas no quinto parágrafo e faz-se a opção 
daquele que melhores resultados fornece para o estudo 
em questão. 

Já na análise do sinal falado, por intermédio do 
«Cepstrum», estuda-se a técnica de «Short-cepstrum», 
no sexto parágrafo, apresentando-se no sétimo os resul- 
tados da aplicação do «Short-cepstrum» a um sinal fa- 
lado, bem como a análise desses resultados. 

No oitavo parágrafo apresenta-se uma aplicação 
importante do estudo do sinal falado feito nos capítulos 


anteriores, utilizando o «Cepstrumi — os vocoders» — 
mostrando-se também diversos tipos possiveis de «vo- 
coders». 

No nono parágrafo mostraremos que a técnica de 
«Cepstrum» deixa de ser útil para sinais falados imersos 
em ruído com pequenas relações sinal-ruído, sendo então 
necessário empregar outras técnicas para a determinação 
do «pitch». A estas técnicas se faz referência no dé- 
cimo parágrafo. 

Chama-se a atenção para o facto de O nosso traba- 
lho ser um estudo experimental, fazendo-se, no entanto 
sempre que possivel, a sua fundamentação teórica. 

Por último, queremos também salientar que a téc- 
nica de «Cepstrum» cujo estudo se faz é igualmente uti- 
lizada noutros dominios, tais como sismologia, radar e 
sonar o de um modo geral todos os sinais que apresen- 
tam ecos (ver o trabalho referido acima, [6]) e ainda 
no tratamento de sinais biológicos, processamento foto- 
gráfico, detecção de sinais imersos em ruido, etc., donde 
se infere a importância desta técnica. 


2 -— O SINAL FALADO 


Vamos descrever neste parágrafo o mecanismo 
completo de geração dos sons falados, dado que tal se 
torna necessário para o desenvolvimento do nosso 
estudo, 

Os sons da voz são praduzidos pela expuisão do ar 
dos pulmões (que são, por assim dizer, a nossa fonte de 
energia vocal) e que é feita de uma forma continua, de 
modo que por analogia com as fontas continuas de ener- 
gia eléctrica lhe chamamos fonte «DC». Esta energia con- 
tirua é transformada em energia pulsariíie (energia «ÃC») 
por um dos três processos seguintes: 


a) Por acção das cordas vocais, a qual transforma 
o fluxo continuo de ar em pulsações de curta 
duração e com uma ocorrência quasi-periódica. 


b) No interior do tracto vocal criam-se turbulências 
de a: por variações de contracções desse mesmo 
tracto vocal. 


c) Suúbitas libertações de excessos de pressão se- 
guidas de um completo encerramento do tracto 
em alguns locais ao longo da sua extensão. 

O sinal acústico «AC» formado por um destes três 

modos é chamado «função de excitação» da cavidade 
vocal e pode tomar três formas diferentes consoante o 
modo da sua formação. Assim temos respectivamente: 


a) Pulsações quasi-periódicas. 

b) Ruído continuo, 

c) Uma pulsação única. 

A maior parte dos sons são gerados a partir de uma 


única destas funções de excitação, Assim, são exemplos 


— — E 


de sons formados pela forma a), que na literatura anglo- 
-saxónica aparecem com a designação de sons «voiced» 
e que poderemos traduzir por sons vocalizados, os sons 
a, e,i, 0 e u. Para exemplos de sons gerados pelas for- 
mas bj e c), chamados sons «unvoiced» (sons não vocali- 
zados) temos respectivamente s, f, z e p, t. 

Alguns sons são gerados a partir de uma mistura 
de dois destes modos e muito poucos utilizam os três 
conjuntamente. 

Para produzir os diferentes sons falados, para cada 
modo de excitação, a função de excitação é filtrada pelo 
tracto vocal, cuja função de transferência (evidente- 
mente no domínio da frequência) depende da posição 
da lingua, dos lábios e outras articulações. No discurso 
esta função de transferência vai, portanto, variando no 
tempo, mas é sempre caracterizada por um certo nú- 
mero de ressonâncias, a que se dá o nome de formantes, 
frequentemente abaixo dos 3000 Hz (valores típicos: 
500 Hz, 1500 Hz e 2500 Hz), como veremos adiante 
mas pode ser já percebido se analisarmos as figuras 4 2 
a 4.4. 

Esquematicamente, todo o mecanismo descrito atrás 
pode ser representado como se segue: 


.—-- — — — — = 


h(t) 


pq — — a 


Re Per )=pltjen(e) 
convergor do "DE" em "AgM, 


FONTE VOCAL 


onde se representou por pi(t) a função de excitação do 
tracto vocal e por h(t) a resposta do tracto vocal ao 
impulso unitário. Os blocos b) e c) representam os pro- 
cessos de formação da função de excitação atrás des 
critos por b) e c) respectivamente (sons «unvoiceda). 
O trajecto que no esquema parte dos pulmões e passa 
pelas cordas vocais é o que produz os sons «voicedh. 
O sinal falado assume, portanto, a forma 


ft) = pl) ht) (0) (2.1) 


ou seja, a convolução de pit) com hit). 

Para o estudo que se segue a partir de agora « até 
ao parágrafo 7, inclusivé, interessa-nos principalmente 
o mecanismo de formação dos sons vocalizados. Para 
estes sons a função p(t) é, como já vimos, uma sucessão 
de impulsos de curta duração produzidos pelo fluxo de 
ar através das cordas vocais as quais permitem a sua 
passagem a intervalos de tempo mais ou menos regulares 
dando origem a ondas de pressão cuja representação 
gráfica está feita na figura 2.1. 


(1) Matematicamente, sabemos que t' é uma variável diferente de t, o que resulta da definicão de autocorrelação YZ) = 


Os; 
=lim —— | 
T-s00 T " 


variável t, tendo sempra presente este facto. 


fl). flt+ 2). dt, donde t = 7, No entanto, por comodidade de exposição, passaremos sempre a usar a 
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ed Ce dO Foda ai Ses LEA ana 


Fig. 2.1 


Estes impulsos de pressão originados nas cordas 
vocais vão excitar os modos de oscilação possíveis no 
tracto vocal que actua assim como uma cavidade de res- 
sonância. 

Note-se que tanto o período T do sinal p(t) pro- 
duzido nas cordas vocais, como a forma do tracto vocal 
e portanto a resposta ao impulso unitário h(t) (ou seja 
os seus modos de ressonância) são variáveis no tempo. 
Pode pois considerar-se o tracto vocal como um filtro 
linear, cuja função de transferência é variável no tempo 
e cuja entrada é uma função «periódica» com um «pe- 
ríodo» igualemnte variável no tempo. É possível, porém, 
admitir hipóteses simplificativas, dado que a função de 
transferência do tracto vocal H(w) = TF[h(t)] (transfor- 
mada de Fourier de h(t)) não varia muito rapidamente 
com o tempo durante alguns «períodos» de p(t) e dado 
que o «período» de p(t) varia lentamente com o tempo, 
de modo que p(t) é uma função quasi-periódica. H(«w) 
e T podem, pois, considerar-se constantes em intervalos 
de tempo contendo vários «períodos» de p(t). Veremos 
em nota no parágrafo 6 a justificação e validade destas 
afirmações. 

Do que se disse sobre a formação de sons «voiced» 
so infere que o sinal f(t) é quasi-periódico (de período 
aproximadamente constante para intervalos de tempo 
n.T, com n inteiro e pequeno) e resulta da convolução 
de p(t) com h(t). Na figura 2.2 apresenta-se um exemplo 
de som «voiced» captado através de um microfone de 
condensador de alta qualidade, filtrado entre O e 2500 Hz 
e amostrado com uma frequência de amostragem de 
5000 Hz, conforme o teorema da amostragem. O conver- 
sor analógico-digital utilizado para fazer a amostragem 
faz parte de um analisador de Fourier HP 54524A, unidade 
periférica de um computador digital HP 2116C que 
efectuou o registo. Desse registo foi retirada uma amostra 
de 128 pontos, ou seja, com uma duração de 25,6 mseg. 
Para o traçado do troço escolhido do sinal registado, 
através de um «plotter», o computador executou previa- 
mente uma interpolação entre os pontos registados (in- 
terpolação linear). Este processo foi utilizado sempre que 
se traçaram gráficos e deve ser sempre tomado em conta 
sem o que se poderão fazer deduções erradas, 


(ajue(ajd = (2)3 


Se for P(«w) a transformada de Fourier do sinal ptt) 
da fonte vocal, 
Plo) = TF[p(t] (2.2) 
e se H(w) for a transformada de Fourier da resposta do 
tracto vocal ao impulso unitário h(t) (função de transfe- 
rência do tracto vocal) 
H(w) = TF[h(t)] (2.3) 
admitindo que ambas as transformadas de Fourier exis- 
tem, teremos: 
F(w) = P(w).H(w).27 (') (2.4) 
sendo 
F(w) = TF[f(t)] (2.5) 
Na frequência, o sinal cujo espectro de frequências 
é P(c), é filtrado pelo tracto vocal. De notar que os 
modos de oscilação do tracto vocal aparecem como picos 
na função de transferência H(«w), chamando-se aqueles 
modos formantes da voz. 
De (2.4) temos 
20 log F(w) = 20 log P(w) + 20 log H(cw) (2.6) 
Analisando agora a figura 4.2 onde se representa 
| F(w) | em dB, verifica-se que se pode explicitar a es- 
trutura formante (a traço-ponto), que pode ser obtida 
de uma maneira aproximada, com bastante rigor, tra- 
cando a envolvente da curva | F(«w) | (curva a traço inter- 


rompido), facto este que nos vai ser muito útil no pará- 
grafo 8. 


5 — NECESSIDADE E UTILIDADE DO «CEPSTRUM» 


Como veremos mais detalhadamente nos próximos 
parágrafos, o problema fundamental que se nos põe é o 
da determinação do «pitch», ou seja, do período T do 
sinal da fonte p(t). É precisamente de um dos processos 
de determinação deste «pitch» que nos vamos ocupar. 

Para determinar esta frequência fundamental de um 
sinal periódico, o método que se seguiria «a priori» seria 
o seguinte: 

Se o período do sinal f(t) é T segundos, o seu es- 
pectro de potência | F(«w) |* é formado por picos espaçados 
de 1/T Hz, correspondentes à frequência fundamental e 
todas as suas harmónicas, o que quer dizer que o es- 
pectro de potência de um sinal falado quasi-periódico é 
«periódico» ao longo do eixo das frequências com «pe- 
riodo» igual a 1/T Hz, ou seja, o inverso do período de 
p(t). 

O método que então utilizaríamos para determinar 
a periocidade do espectro de potência, 1/T Hz, seria 
calcular a transformada inversa de Fourier do espectro 
de potência. Esta transformada teria um pico correspon- 
dente a esse «período». 


(2) A partir de agora não voltaremos a apresentar o factor 2, por não ser relevante, 
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Este espectro do espectro de potência é conhecido 
pelo nome de função de autocorrelação do sinal no 
tempo f(t) — Teorema de Wiener-Khintchine. 

Sendo Ye tc) a função de autocorrelação de fit) 
teremos, por definição 


2 q 
ve (o) = — rea) Leto 2 | (3.1) 
$ 
8 
; “gh 1 
rlo=, Polo —tustt+ola (3.2) 
po 00 Jotta 


Como o espectro de f(t) é 
F(w) = P(w).H(c) (3.3) 
teremos consequentemente 


TE" [IFlo)l2]=TES [1P(u) 2] *TE- 
[HH(e) 1*] (3.4) 


onde * representa a função convolução. 


Será então: 

ta MO) = Yo 1)" Pg te) (3.5) 
com 

CS TEA [1 P(w)1*] (3.6) 

yn (2) = TES [UH(o) 17] (3.7) 


Assim chegémos à conclusão de que a autocorrela- 
ção de f(t) é a convolução das autocorrelações de pít) 
e de h(t). Portanto os efeitos do som da fonte e do 
tracto vocal são convolucionados um com o outro na 
função de autocorrelação. 

Deste facto resultam várias dificuldades nara a de- 
terminação do «pitch», entre os quais temos: 

Para certos sinais falados a função de transferência 
de tracto vocal é tal que a autocorrelação Y; (7) tem 
um grande número de picos. Noutros casos os picos 
apresentam uma grande largura ou ainda a autocorrela- 
ção apresenta um patamar. 

Em qualquer destes casos a determinação do «pitch» 
é difícil, senão mesmo impossível, quer por dificuldade 
em detectá-lo, quer em defini-lo com precisão. 

Consideremos por exemplo o sinal falado da figura 
3.1 que é o sinal que se mostrou atrás na figura 2.2. 
A sua autocorrelação não nos dá uma informação precisa 
sobre o «pitch», como se vê na figura 3.2, vor o pico 
ser muito largo, apresentando, portanto uma grande im- 
precisão. 


e KIM SM 


Fig. 3.1 —- Sinal falado filtrado entre O e 2,5 kHz (frequância de 
amostragem: 5 kHz) Número de pontos registados: 128; 
T = 8,4 mseg 
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Fig. 3.2 —- Autocorrelação do sinal falado representado na fig. 
3.1; máximo para t = 8,6 mseg 
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Fig. 3.3 — Sinal falado, filtrado entre O e 2,500 Hz (frequência 
de amostragem: 5 kHz) Número de pontos renistados: 512; 
T=7 mseg 
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Fig. 3.4 — Autocorrelação do sinal da fig. 3.3 


Para o sinal da figura 3.3, cuja frequência funda- 
menta! é facilmente calculada por exame directo, a auto- 
correlação é a apresentada na figura 3.4 e tem, como 
se vô uma multiplicidade de picos. Neste caso não seria 
difícil obter informações sobre o período T, se os valores 
do dois picos consecutivos não fossem, frequentemente, 
próximos um do outro. Os picos poderão mesmo ter 
ocorrência não periódica, dificultando a determinacão 
do «pitch». Mas mesmo que a ocorrência seja periódico, 
tomando em conta que por vezes se dá a duplicação do 
periodo T, o que acontece de uma maneira progressiva 
numa mesma amostra de sinal falado poderemos ter dois 
«pitch», um com o valor duplo do outro, situação esta 
quo importa igualmente detectar. Com a autocorrelação 
apresentando picos múltiplos periódicos não sabemos 
nunca se estamos na situação descrita ou não, tornando 
a análise complicada. 

Sabemos, por outro lado, que para sinais de banda 
estreita a autocorrelação não fornece resultados satisfa- 
tórios [6]. Para estes sinais a utilização do «tapering» 
seguida do «cepstrum», como se explica no parágrafo 6, 
resolve o problema por razões que então veremos. No 
entanto não conseguimos descobrir nenhum sinal falado 
que fosse de banda estreita (embora o pudéssemos obter 
artificialmente por filtragem, achámos incorrecto fazê-lo, 
porque corresponderia a uma situação em que se perderia 
a inteligibilidade da fala), pelo que não apresentamos 
qualquer exemplo deste último caso. 

Problemas de vária ordem, como por exemplo, 
pequenas relações sinal-ruído, distorções não lineares, 
mudançasde fase, etc., tornam o problema da determina- 
ção do «pitch» de um sinal cujas características estatis- 
ticas não são conhecidas «a priori» num problema com- 
plicado [2]. Nos parágrafos 9 e 10 faremos uma refe- 
rência resumida a estes casos. 
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Do que ficou dito (excepto no narágrafo atrás que 
exige outros métodos) se infere da necessidade e utili- 
dado do «cepstrum», o qual vamos passar a estudar em 
pormenor no próximo parágrafo. 


4 — «CEPSTRUM» 


Vimos, portanto, que a determinação do «pitch» 
através da autocorrelação não é em geral satisfatória. 


Considere-se, então, que antes da determinação da 
periodicidade 1/T Hz do espectro de potência de f(t), 
| F(w) |* se emprega a função logaritmo. 


| F(o) | =1P(o) FF. 1H(o) [É (4.1) 


log | Flw) 1* = logl P(w)|* + log | H(e) |“ (4.2) 

Vemos que deste modo os efeitos do tracto vocal 
e da fonte vocal so adicionam (a transformada inversa de 
Fourier do uma soma é igual à soma das transformadas 
inversas de Fourier das parcelas) em vez de se convolu- 
cionarem, como acontecia não empregando a operação 
intermédia de logaritmo. 


Note-se que não teremos de trabalhar necessaria- 
mente com espectros de potência. Se utilizarmos os es- 
pectros de frequência teremos 


[F(o)L=1P(o)1.!H(0)| (4.3) 


log | Flw) | = log l P(w)|l + log | H(w) | (4.4) 
e O tratamento posterior é análogo ao que vamos desen- 
volver usando espectros de potência, pois 

log | P(w) |" = 2log | P(w) | (4.5) 
e portanto há apenas uma mudança de escala na passa- 
gem de um processo para outro. Como veremos no pa- 
rágrafo seguinte, é com este último processo que traba- 
lharemos, pois, para os meios de cálculo que utilizámos 
conduz-nos a erros menores. 


A importância de os efeitos do tracto vocal e da 
fonte vocal serem aditivos pode ser facilmente visionado 
nos exemplos seguintes. 


Na figura 4.1 está representado o módulo do es- 
pectro de frequências do sinal da figura 3.3 e na figura 
4.2 o seu logaritmo (ver final do parágrafo 2 nara a in- 
terpretar). Nas figuras 4.3 e 4.4 temos respectivamente 
o espectro de potência do sinal da figura 3.3 e o seu 
logaritmo. Tanto no logaritmo do módulo de F(w) como 
no logaritmo de | F(w) |” o efeito do tracto vocal é pro- 
duzir uma oscilação de baixa «quefrência» (anagrama da 
palavra frequência, sugerido pela troca de domínios: 
quefrôncia = frequência de uma função com domínio na 
frequência [5]), enquanto que a periodicidade da fonte 
vocal vai dar origem a uma «periodicidade» de 1/T Hz, 
ou seja, uma oscilação de alta frequência. 
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Para dete:-minar esta periodicidade 1/T Hz podemos 
utilizar um dos métodos próprios e bem conhecidos para 
o fazer, como sejam o espectro de frequências ou o es- 
pectro de potência. 

Mais rigorosamente vamos determinar o espectro 
de quefrências (ou o quadrado do seu módulo) da 
função 


log | Flw) 1|* = log I P(w) |* + log | H(w) |* (4.6) 
ou da função 
log | Flw) 1 = loglP(w) | + log Ht(w) |. (4.7) 


O espectro de quafrências, ou seja, a transformada 
inversa de Fourier de qualquer das funções (4.6) ou 
(4.7) apresenta um pico agudo correspondente à alta 
quefrência, o «pitch», e picos diversos para valores 
baixos da quefrência que correspondem à estrutura for- 
mante. É o que se vô na figura 4.5, onde está represen- 
tada a transformada inversa de Fourier do logaritmo do 
módulo do espectro de frequências do sinal da figura 
3.3 e na figura 4.6 onde se representou a transformada 
inversa de Fourier do logaritmo do espectro de potência 
do mesmo sinal. 

Note-se que sendo 


log | F(w) |* = 2.1og | F(w) 


TF (loglF(w) I* ) = 2.TF“ | logl F(w) 1] (4.8) 


donde as funções são as mesmas, aparte o factor multi 
plicativo 2. 


Fig. 4.1 —- Especiro de frequência (módulo) do sinal da fig. 3.3, 
1 


T = 143 Hz 


estrutura firma 
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Fig. 4.2 —- logaritmo do módulo do espectro de frequências do 
sinal da fig. 3.3. (ultimo ponto sem significado) 
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Fig. 4.3 — Espectro de potência do sinal da fig. 3.3 
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Fig, 4.4 —log | Fl,) |! do sinal da fig. 3.3 
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Fig. 4.5 — «Cepstrum» de frequência do sinal da fig. 3.3 (TF -' 
"log | Ft) 13) Pico máximo para = 7 mseg (excluindo os 
picos em torno da origem). 


Fig. 4.6 —- «aCepstrum» de potência do sinal da fig. 3.3 (TF -' 
"log (1 Fls) E) 1) Pico máximo: Z = 7 mseg (excepto picos em 
torno da origem) 


Façamos notar ainda que este método, além de ex- 
tremamente útil para a determinação do «pitch», nos per- 
mite também conhecer as quefrências das formantes (que 
não são mais que os seus periodos) e é aliás o que for- 
nece melhores resultados. E isto é muito útil para a sín- 
tese da voz como se faz nos «vocoders», como veremos 
no parágrafo 8. 

Resta acrescentar que o ponto na origem (7 = O) 
não é relevante para os nossos fins (determinação das 
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formantes e, principalmente, do «pitch), pois, como é 
óbvio, corresponde à componente contínua de log | Flw) | 
(«Cepstrum» de frequência) ou de log | F(w) |” («Ceps- 
trum» de potência). 


E — «CEPSTRUM» DE POTÊNCIA E «CEPSTRUM» DE 
FREQUÊNCIAS 


Para tornar os resultados mais evidentes podem 
elevar-se ao quadrado os módulos das transformadas in- 
versas de Fourier das figuras 4.5 e 4.6. É ao quadrado 
do módulo da transformada inversa de Fourier do loga- 
ritmo do espectro de potência 


ITE flogl Flo) IF) 


que Tukey chama «CEPSTRUM» [3], que é um anagrama 
da palavra espectro resultante da troca de domínios 
atrás referido. No entanto nós chamaremos «cepstrum» 
tanto a esta função como a TF, (logl F(w) 1), bem 
como a TF-* (logl F(w) |"). O que nós vamos utilizar, 
contudo, será o segundo que chamaremos de «cepstrum» 
de frequência para o distinguir, quer do «cepstrum» de 
Tukey, quer do terceiro, a que chamaremos «cepstrum» 
de potência. No entanto, como vimos atrás, eles estão 
intimamente relacionados, dado que este último não é 
mais que duas vezes o «cepstrum» de frequência e o 
«cepstrum» de Tukey é quatro vezes o quadrado do 
«cepstrum» de frequência, o que para os nossos fins 
tanto faz. Com efeito: 


«Cepstrum» (de Tukey) = 

= [TF-* [loglF(uw)1*)])* 
[TF-* | 2loglF(wo)1)]: 

= 12.TF-* (logl Flw) lj1* 

= 41ITF-* iloglF(wo)1)]' 

= 4 | «cepstrum» de frequência |" 


De notar que a quefrência se exprime em ciclos/ 
/(ciclos/seg) ou, simplesmente, em segundos. 

Analisemos, por exemplo, os resultados obtidos 
com o sinal da figura 3.3 obtidos no capítulo anterior. 

O sinal representado tem a duração de 102,4 mseg 
e apresenta quinze picos nesse intervalo de temno pelo 
que a sua frequência fundamental é de cerca de 150 Hz 
ou seja T = 7 mseg. O espectro de frequências, bem 
como o logaritmo do seu módulo, têm, portanto uma 
periodicidade de 1/T Hz = 150 Hz, o que se pode con- 
firmar nas figuras 4.1, 4.2, 4.3 e 4.4. Nas figuras 4.5 e 
4.6 vemos que os picos de amplitudes máximas, numa 
e noutra (exceptuando, evidentemente os picos em torno 
da origem que, como vimos, correspondem à estrutura 
formante) se verificam para Z = 7 mseg. 

Analizando mais atentamente as figuras referidas 
notamos que existem outros picos de amplitudes mais 
reduzidas para 7 = 14 mseg e Z = 21 mseg. Se nos 
recordarmos do modo como foram obtidas, vemos que 
o último passo dado foi o cálculo de uma transformada 
(inversa) de Fourier, que além de determinar a «fre- 
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quência» fundamental do logaritmo do módulo da trans- 
formada de Fourier de f(t) (fig. 4.5) ou do logaritmo do 
quadrado desse módulo (fig. 4.6). determina todas as 
suas «harmónicas». Em geral os picos correspondentes 
a estas «harmónicas» estão dissimulados na irregulari- 
dade do cepstrum e não nos trazem, por isso quaisquer 
problemas. 

Vejamos agora as razões que nos levaram a utilizar 
o «cepstrum» de frequência e não o «cepstrum» de po- 
tência. Para chegar a este último calculamos dois qua- 
drados de módulos de espectros, o que não acontece 
na determinação do primeiro. Analisando, por exemplo, 
r módulo do espectro representado na figura 4.1, vemos 
que ele é constituído por um conjunto de picos de valor 
elevado numa dada faixa de frequências (entre zero e 
1,3 kHz, aproximadamente) e que fora desta faixa e entre 
os picos tem valores baixos mas que nalgumas zonas de 
frequências têm alguma importância. Ao calcularmos o 
quadrado desta função reparámos que o computador 
cometia muitos erros resultantes das quantificações que 
um computador digital necessariamente faz, erros esses 
tanto mais importantes quanto menores forem os valores 
em jogo. 

Assim, vemos na figura 4.3, o pico de valor mais 
elevado é posto mais em evidência, quando se calcula 
o quadrado, como seria de esperar, mas os quadrados 
dos valores baixos atrás referidos são arredondados para 
zero por se tornarem demasiado pequenos. O espectro de 
potência passa assim a ser, erradamente, um espectro 
de bandas, enquanto o espectro de frequência não o era. 
Este facto conduz a erros importantes quando se calcula 
o «cepstrum» de potência. 

Analisemos, por exemplo, as figuras 5.1 e 5.2. Na 
primeira está traçado o «cepstrum» de frequência da 
secção do sinal falado «ATA» da figura 6.1, entre t = 
= 76,8 mseg e t = 102,4 mseg. Essa secção não foi ob- 
tida multiplicando o referido sinal pela função rectan- 
gula: centrada em t = 89,6 mseg (ponto médio do inter- 
valo considerado), mas sim pela janela representada na 
figura 6.2. A razão desta escolha será explicada no capí- 
tulo seguinte, não interessando neste momento, pois o 
que nos propomos é mostrar que o «cepstrum» de potên- 
cia de um sinal qualquer pode conduzir a erros impor- 
tantes utilizando os meios de cálculo que nós usámos. 
Na segunda está representado a TF-! ( log | F(w) 17), 
que é o nosso «cepstrum» de potência. Já vimos atrás 
que este é o quádruplo do quadrado do módulo daquele, 
o que, é fácil de ver, não se verifica. O «cepstrum» de 
potência, devido às razões atrás apontadas, não conduz 
a resultados satisfatórios. 


E imseg) 


Fig. 5.1 — «Cepstrum» frequência da secção de sinal represen- 
tado na fig. 6.1 entre t= 76,8 mseg e t= 102,4 mseg 


mostrará a sua validade. 
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Fg. 5.2 — «Cepstrum» de potência da secção de sinal referida 
na figura anterior 


6 — «SHORT-CEPSTRUM»; ANÁLISE TEÓRICA. 
«TAPERING» — JANELA NO TEMPO 


A técnica de «Cepstrum» pode aplicar-se a um sinal 
de qualquer comprimento no tempo e o que se fez até 
agora foi determinar o «cepstrum» de um sinal com uma 
determinada duração no tempo, mas tomado sempre 
como um sinal inteiro. 

Este é o método utilizado no estudo de sinais sís- 
micos, e de um modo geral, de todos aqueles em que o 
que se pretende determinar é a existência de um ou vá- 
rios ecos e os respectivos tempos de ocorrência, relati- 
vamente ao sinal «principal». 

Mas no caso do sinal falado, a periodicidade do 
sinal emitido pela fonte vai variando no tempo, e há 
mesmo porções do sinal falado que são aperiódicas. 

Considere-se, por exemplo, o sinal falado da fi- 
gura 6.1, que foi recolhido, amostrado, registado e tra- 
çado nas mesmas condições dos sinais das figuras 3.1 
e 3.3. Trata-se da palavra «ATA» (do verbo atar) dito 
por uma pessoa do sexo masculino, com a duração apro- 
ximada de 200 mseg (o filtro utilizado para filtrar tanto 
este sinal como os dois referidos atrás foi um filtro B.W. 
passa-baixo O — 2500 Hz). Este sinal é constituído por 
três partes distintas: a primeira (O a 50 mseg, aproxima- 
damente), diz respeito a um som quasi-periódico («voi- 
ced») a que corresponde um «pitch», eventualmente va- 
riando no tempo; a segunda (aproximadamente entre 5O 
e 150 mseg), é um som não periódico («unvoiced») e 
corresponde à pronunciação do som T, durante o qual 
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Fig. 6.1 —- Sinal falado «ATA» dito por um indivíduo do sexo 
masculino, filtrado entre O e 2,5 kHz (frequência de amostragem: 
5 kHz). 


não há admissão de ar através das cordas vocais e ao 
qual não há pois que atribuir qualquer «pitch»; a terceira 
parte é idêntica à primeira. 

Como a característica «voiced-unvoiced» do sinal 


falado vai variando no tempo, bem como o período do 
sinal da fonte, convém-nos analizar o sinal por secções, 
isto é, considerando primeiramente um troço do sinal 
entre t= 0 et=2t' e obter o «cepstrum» dessa porção 
de sinal, analizindo se ela é a«voiced» ou unvoiced» e no 
primeiro caso qual o «pitch» que lhe corresponde, Em 
seguida consideramos outra porção, que no nosso caso 
do sinal da figura 6.1, é entre t= t' e t= 3t' repetindo 
a análise feita para a primeira, e assim sucessivamente. 
De notar que a análise passa a ser feita por saltos de 
At, no nosso caso At =t”, e que as diversas porções 
(todas elas com a mesma duração — no nosso caso 2t') 
se podem sobrepor. Esquematicamente (para o nosso 
caso): 


|.* porção 3.“ porção 5.“ porção 
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É neste ponto que é necessária a hipótese feita na 
página 8, de que tanto o «pitch» do sinal falado, como 
a função de transferência do tracto vocal H(«w) variam 
lentamente com o tempo. Com efeito a divisão do sinal 
falado em porções, como atrás é referido, destina-se a 
atribuir a cada uma delas a classificação de «voicedy 
ou «unvoiced» e. no primeiro caso, um determinado 
«pitch», que se considera constante em cada uma delas. 
Ora, se as variações do «pitch» e da função de transfe- 
rência H(w) com o tempo fossem muito rápidas, isto 
seria manifestamente impossível de fazer. Não seria 
possível atribuir a uma dada porção do sinal a caracte- 
rística de «voiced» ou «unvoicedys e muito menos um 
«pitch» determinado. 

A duração das porções (2t') bem como os saltos 
A t deevm ser escolhidos conforme as velocidades de 
variação do «pitch» e de H(w). 

No nosso caso, foram julgados suficientes saltos e 
porções de 12,8 e 25,6 mseg respectivamente, que, como 
veremos no próximo parágrafo conduziram a resultados 
aceitáveis. Aí veremos igualmente que é válida a hipó- 
tese das pequenas variações para o «pitch» e para H(«).(*) 

Ão seccionarmos o sinal falado, o que estamos real- 
mente a fazer é a multiplicá-lo por uma função rectangu- 
lar. Este facto vai dificultar-nos os cálculos pois vai-nos 
introduzir na porção de sinal a analisar frequências não 
existentes no sinal inicial. Vai haver um enriquecimento 
especialmente no que diz respeito às altas frequências. 

Vamos analisar, seguidamente, qual é este efeito, 
como se dá e como o poderemos resolver utilizando, 
em vez de uma função rectangular, uma outra janela 
no tempo especialmente escolhida. 

Sendo o sinal a analisar 


f(t) = plt) * htt) (6.1) 


mostrará a sua validade, 


multipliquemo-lo por uma dada função wít): 
g(t) = Hit.wlt) = (plt) * h(t)) wit) (6.2) 


wit) é a função utilizada como janela, que tanto pode 
ser a função rectangular como outra qualquer. 

O sinal p(t), da fonte, representado graficamente 
na figura 2.1 pode ser expresso por uma soma de uma 
sério infinita do sinais idênticos po(t), repetidos cada 
T segundos, sendo T o «pitch» de pít), isto é: 


+" 00 
p(t) = = polt-nT) (6.3) 
n=-2 
+= 00 
=polt)* E s (tnT) (6.4) 
n=-0€ 


Substituindo em (6.2) teremos: 
+ 00 
att) = [lpolt)* E o s(tnT)) * h(t))wit) (6.5) 


n=-o 


e a transformada inversa de Fourier deste sinal é: 


+ 27 
G (0) -| Polo). E leon) to | * ww) (6.6) 
En A 
n=-c 
+90 27 
= | H (6). Polo). sen) * ww) (6.7) 
Er E 
n=-00 
+90 2% 2% 2a 
pê sfMto Deum ” ). S(w-n z pp E wlw) (6.8) 
n=-00 


Deste modo, o espectro de frequências do sinal fa- 
lado é formado por impulsos espaçados de 1/T Hz e com 
amplitude H(n 27/T). Pon 27/T) se o sinal falado não 
foi limitado no tempo (w(t) “- O, excepto para valores 
isolados do t), o que aliás já tinhamos visto atrás. 

Não sendo o sinal limitado no tempo é limitado na 
frequência, ou seja, O seu espectro de frequências é nulo 
para valores de «w tais que | w | > — a 

Sendo, o sinal falado, limitado no tempo pela ja- 
nela w(t), o efeito desta janela, traduz-se, no espectro 
de frequências, pela convolução do espectro de frequên- 
cias do sinal como se ele fosse ilimitado no tempo e por- 
tanto, limitado na frequência, com o espectro de fre- 
quências de w(t), ilimitado na frequência. Quer isto dizer 
que o espectro de g(t) = f(t).w(t) não é limitado na fre- 
quência (vão-nos aparecer frequências estranhas a f(t) 
devido à janela no tempo w(t) ) sendo ainda os picos 
que apareciam em F(«w) alargados e tornados com a 
forma de Wo). 

O espectro de frequências da janela no tempo, fun- 
ção rectangular é: 


(3) Aqui fica mostrada a necessidade destas hipóteses como ss havia prometido no parágrafo 2. No próximo parágrafo se 
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Ty 
sen w-—— 
2 


= 
Wiw) =—É (6.9) 
2 w 


2 


sendo Ty o intervalo de tempo em que a janela toma 


+ fmsegt 
-T 8 Te 


mt 
Fig. 6.2 — Janela no tempo wlt) = 0,54 + 0,46 cos re) 
w 


ti < Tw 


wit) = O Lt1>s Tw 


valores não nulos, o que vem ilustrar o que foi dito 
atrás. 


Assim, se nos limitarmos a seccionar o sinal, o que 
equivale a multiplicá-lo pela janela no tempo — função 
rectangular — os picos de F(«w) são alargados, tomando 
a forma da função cuja expressão é (6.9), do que resulta 
uma distorção nos resultados, impedindo, muitas vezes 
a determinação do «pitch». Este é o rpoblema que se põe 
para os sinais de banda estreita, como já referimos atrás, 
pois é para estes que, naturalmente, se verificam as 
maiores distorções no seu espectro de frequências e, 
portanto, no seu «cepstrum», não permitindo, geralmente 
que se detecte o «pitch» de tais sinais. 


O problema que se nos põe agora é escolher conve- 
nientemente uma janela no tempo w(t), de modo que, 
não sendo o espectro de frequências de g(t) limitado na 
frequência, possa contudo ser aproximadamente consi- 
derado como tal. Se o espectro W(«w) da janela for apro- 
ximadamente limitado na frequência, os picos de F(cw) 
não serão muito arredondados, não haverá um excessivo 
enriquecimento em frequências que não pertençam a fít) 
e obtêm-se os resultados desejados. 


Por isso a janela no tempo deve ter um espectro 
de frequências com um lobo principal o mais largo pos- 
sível e com os lobos secundários com a menor amplitude 
possível. Estas duas exigências são mutuamente exclu- 
sivas, pelo que é necessário tomar uma solução de com- 
promisso. 


Uma boa janela no tempo, tomando em atenção 
o que se disse anteriormente, é a que tem por expressão 


qt 
wit) = 0,54 + 0,46 cos(—) 
Ty 


para ti ST 


wíit) = 0 para IT 


em que 2 Tw é 0 intervalo de tempo da secção do sinal 
e analizar, e que tem a forma representada na figura 6.2. 
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Fig. 6.3 — Secção entre t=0 e t= 25,6 mseg do sinal falado 
«ATA» da fig. 6.1 
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Fig. 6.4 — Sinal da fig. 6.3 multiplicado pela janela da fig. 6.2 
(Fw = 12,8 mseg) (Abreviadamente: com «tapering>) 
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Fig. 6.5 — Espectro de frequências da janela da fig. 6.2 (módulo 
em dB) 


Nas figuras 6.3 e 6.4 estão representados respecti- 
vamente a secção entre O e 25,6 mseg do sinal da figura 
6.1 e o produto desta pela janela representada na fi- 
gura 6.2. 

A janela descrita atrás apresenta uma relação lobo 
principal-lobos secundários de 25 dB, como se pode ver 
na figura 6.5, onde está traçado o espectro de frequên- 
cias dessa janela. 

Resta acrescentar que com esta janela se obtêm 
bons resultados como se verá no parágrafo seguinte, 
mas que é possível chegar a bons resultados igualmente 
utilizando outras janelas, como acontece por exemplo 


em [61. 


7 — «SHORT-CEPSTRUM» DE UM SINAL FALADO: 
ANÁLISE DOS RESULTADOS 


Vamos apresentar, neste capítulo, os resultados ob- 
tidos pela aplicação da técnica atrás descrita, e que é 
conhecida pela técnica de short-cepstrum, ao sinal falado 
«ATA» da figura 6.1. A janela utilizada, e que foi definida 
no parágrafo anterior, tem uma duração de 25,6 mseg 
como então se disse, e o «short-cepstrum» foi feito 
por intervalos de 12,8 mseg. Dado que o sinal completo 
tem uma duração de 204,8 mseg, obtivemos 15 «ceps- 
tra». De referir, antes de procedermos à sua análise, que 
o sinal foi dito por uma pessoa do sexo masculino, cap- 
tado, filtrado, amostrado e registado nas mesmas condi- 
ções que o sinal da figura 2.2 (captado por um micro- 
fone de condensador de alta qualidade, filtrado por um 


filtro B.W. entre O e 2500 Hz, amostrado pelo conversor 
AD do Fourier Analyser, unidade periférica do computa- 
dor digital HP 2116C que registou o sinal e, em seguida 
executou todos os cálculos). 


Nas figuras 7.1 a 7.9, inclusive, apresentam-se três 
secções de sinal, os seus espectros logarítmicos e os 
seus «cepstra» de frequência, que a partir deste momento 
passaremos a chamar simplesmente de «cepstra» dado 
que o que utilizamos sempre são os de frequência. 


As duas secções das figuras 7.1 e 7.4 são as duas 
primeiras, respectivamente entre t= 0 e t= 25,6 mseg 
e entre t= 12,8 mseg e t = 38,4 mseg e que, como já 
dissemos se sobrepõem. Tanto numa como noutra é pos- 
sível, por análise directa, determinar o «pitch», pelo 
menos de uma forma aproximada, e ter uma ideia da sua 
estrutura formante por análise do espectro logarítmico. 
Nos «cepstra» vemos imediatamente que o «pitch», tanto 
num como noutro é de 7,4 mseg e que, portanto, não 
variou. À secção de sinal da figura 7.7 é de uma região 
do sinal em que este é «unvoiced», como é evidente, não 
apresentando, por isso, qualquer «pitch», como se torna 
explícito por análise do seu «cepstrump». 


) an R| À 


Fig. 7.1 — Secção entre t=0 e t= 25,6 mseg do sinal falado 
«ATA» dt fig. 6.1 


Fig. 7.2 —- Módulo (em dB) do espectro do sinal da fig. 7.1 


com «tapering+» 
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Fig, 7.3 —- «Cepstrum» do sinal a que se refere a figura acima 
c/«tapering») Pico para Z = 7.4 mseg 
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Fig. 7.4 — Sinal falado «ATA» entre 1=12,8 mseg e 1=38,4 mseg 
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da 


Fig. 7.5 — Módulo (em dB) do espectro do sinal da fig. 7.4, 
c/«tapering» 


Fig. 7.6 — «Cepstrum» do sinal a que se refere a figura anterior 
(c/«tapering») (Pico cortado por ser irrelevante) 
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Fig. 7.8 —- Módulo (em dB) do espectro de frequências do sinal 
da fig. 7.7 c/«tapering» 
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Fig. 7.9 — «Cepstrum» do sinal a que se refere a figura anterior 
c/tapering 


Na figura 7.10 estão representadas todas as sec- 
ções consideradas e que se apresentam para uma crítica 
dos resultados obtidos. Na figura 7.11 estão representa- 
dos o espectros logarítmicos dessas secções do sinal, 
bem como o «short-cepstrum». Resta-nos fazer uma aná- 
lise dos resultados obtidos. 

Os quatro primeiros «cepstra» apresentam picos 
bem nítidos para Z = 7,4 mseg correspondentes ao 
«pitch» T = 7,4 mseg. No terceiro «cepstrum» os seus 
valores para 2= 7,2 e €= 7,4 mseg são muito próximos. 
Como o que está traçado entre eles não tem significado, 
pois que, como o computador é digital, os resultados são 
apresentados sob a forma digital e o que está traçado 
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entre eles é uma interpolação linear feita para o «plotter» eita 
poder traçar continuamente o gráfico, o que temos efec- y ass Eri 
tivamente são dois pontos de amplitudes semelhantes. o a 
Podemos considerá-los como fazendo parte de um pico 
largo ou como sendo dois picos separados, mas em 
ambos os casos a ilacção a tirar é a mesma: dentro do 
intervalo considerado o «pitch» tem uma variação rápida oii é 
e que, portanto, este intervalo entre t= 25,6 e t = 51,2 feios A 
mseg tem, por assim dizer, dois «pitch». 
Note-se, a propósito do que acaba de dizer-se, que 
a janela no tempo atrás considerada, e que foi usada 
neste exemplo, «olha» principalmente para a região cen- 
tral do sinal, dado que as regiões laterais são de pequena 
amplitude. Disto resulta que o pico que aparece no = 
«pitch» nunca é muito largo dado que o «pitch» que é 
predominantemente detectado é o da região central do 
intervalo de 25,6 mseg de que se fez o «cepstrum». 
É possível, no entanto, que por vezes apareçam em torno 
do pico principal, picos de pequena amplitude correspon- DD 
dentes a uma variação do «pitch», ou um pequeno alar- E 
gamento do pico, devido à mesma razão, o que não pa- 
reco ter importância de maior. Se o pico aparecer muito 
alargado, ou se aparecerem vários picos de amplitudes 
semelhantes é peferível fazer uma análise mais fina, uti- 
lizando intervalos mais pequenos e fazendo o «short- nda 


-cepstrum» por saltos menores. se Ei FE ERA 
No quarto intervalo, o cálculo do «cepstrum» con- 

duz-nos a um pico para Z = 7,4 mseg, mas este pico 

é já de pequena amplitude porque, como podemos ver 

na figura 7.10, estamos numa secção em que se faz a 

transição da região em que o som é «voiced» para uma 2 

região em que o som é «unvoiced». 

AO Ah men q 
AN 140,8 


(Za9m)3 
sr or Sh n » Db) 


Jasu)ogóvinp + A+ as 
172,2 
“TS | 
——— qu pp pr duração (nseg) 
[o] ç (0 18 0 13 


Fig. 7.10 —- Sinal falado «ATA» por secções de 25,6 mseg t(maeg) 
sobrepostas Fig. 7.10 — (continuação) 
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Fig. 7.11 —- Espectro logaritmico e «short-cepstrum» do sinal 
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Fig. 7.11 —- (continuação) 
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Fig. 7.11 — (continuação) 
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Nc quinto intervalo estamos já, nitidamente, dentro 
desta região, não aparecendo, consequentemente, nenhum 
pico correspondente a qualquer «pitch». Nesta região os 
picos mais importantes são os que correspondem às fo-- 
manto e que se veem para pequenos valores de z Estas 
formantes são as que modelam, nos sons «unvoiced», 
ruídos produzidos na cavidade vocal. Este assunto é tra- 
tado mais em pormenor no parágrafo seguinte, sendo 
suficiente notar, por ora, que a determinação destas for- 
mantes é tão importante como a des formantes dos sons 
«voiced». 

Nas secções 6 a 10, inclusive, os sinais correspon- 
dem sempre a sons «unvoiced» e, nor isso, os resultados 
são idênticos aos da secção 5, 

Na secção 11, voltamos a estar numa região de 
transição, mas desta vez, de um son «unvoiced» para 
um som «voiced», bem como no intervalo 12, mas em 
ambos os picos são “ão pequenos que se confundem 
com as irregularidades do «cepstrum». 

Nas secções 13 e seguintes, os picos voltam a apa- 
recer de novo nitidamente. É possível notar uma varia- 
ção lenta do «pitch» até ao final do sinal, pois os picos 
nas regiões 13, 14 e 15 correspondem respectivamente 
at = 7,6,€ = 71,8 ez = 8,0 mseg. 

Resta acrescentar que em todos os esnectros loga- 
rítmicos aparece um último ponto (e a interpolação li- 
near entre eles e o penúltimo) que não tem qualquer sig- 
nificado, pois é um erro de cálculo devido ao processo 
de cálculo utilizado pelo Fourier Analyser. 

O processo aqui utilizado pode agora ser usado para 
estudar outro sinal qualquer. 
sentar exemplos com outros sinais, o que não fizemos 
por não dispormos de tempo. 

Seria interessante obter «cepstra» para outros 
sinais e verificar alguns fenómenos conhecidos, como 
por exemplo, o facto de o «pitch» de um dado som 
«voiced» dito por um indivíduo do sexo feminino ser me- 
tado do «pitch» do mesmo som dito por indivíduos do 
sexo masculino, a duplicação do «pitch» nalguns sons, 
como po; exemplo o som inglês «...ing», etc. 


8 — APLICAÇÃO: «VOCODERS» — ANÁLISE ESPEC- 
TRAL E SÍNTESE DA VOZ 


Neste parágrafo vamos ver uma aplicação do estudo 
feito nos parágrafos anteriores. 

Esta aplicação traduz-se no «Vocoder», designação 
derivada das palavras «voice» e «coder». 

Inicialmente o «vocoder» foi concebido como um 
meio para uma eficiente transmisão do sinal falado, es- 
pecialmente para longos e caros circuitos telefónicos. 

Na realidade, numa transmissão vulgar de telefonia, 


a banda ocupada por uma chamada é aproximadamente 
de 3000 Hz, dos 200 aos 3200 Hz. Esta banda é sufi- 
ciente para assegurar uma boa qualidade e inteligibili- 
dado da voz ('). Contudo Dudley mostrou apenas uma 
falados poderiam ser transmitidos utilizando anenas uma 
fracção dessa banda se métodos próprios de código 
fossem usados. O «vocoder» foi assim o primeiro projecto 
para realizar essa economia. 

Actualmente a voz sintética dos «vocoders» tem-nos 
resolvido os problemas da qualidade da voz resultantes 
do elevado nível do ruído, da distorção não linear, da 
distorção de frequência, etc. 

As aplicações dos vocoders têm aumentado e alar- 
gado o seu primitivo âmbito. Actualmente, entre outras 
aplicações, ajudam os surdos a «sentir» os sons da voz 
(tactile-«vocoders»), assim como têm tido larga aplica- 
ção no ensino de línguas estrangeiras. 

Uma das suas mais importantes aplicações, tem sido 
no armazenamento do sinal falado. 

Para o estudo dos «vocoders» é necessário entrar 
em conta, quer com a análise espectral da voz, quer com 
a sua síntese, O que nos leva a estudar as propriedades 
do ouvido humano. 

É necessário considerar as suas três principais pro- 
priedades no que diz respeito à sua sensibilidade ao som 
falado: 


1) O ouvido constitui um analisador espectral em 
«short-time». 


2) Para uma percepção não estereofónica, o ouvido 
é relativamente insensível à fase. 


3) O ouvido é excelentemente sensível ao «pitch». 

No espectro do ouvido ser insensível à fase, deve-se 
precisar que a fase é um factor relativamente sem im- 
portância, tendo alguma influência na qualidade da voz, 
mas sem efeito na inteligibilidade da mesma quando se 
trata de sinais falados não estereofónicos. 

Deve-se notar contudo que distorções de fase cor- 
respondentes a distorções de atrazo superiores a 50 mseg 
já têm efeitos notáveis na qualidade da voz. 

O facto do ouvido ser bastante sensível ao «pitch», 
deve-se ao grande poder descriminatório na frequência 
da percepção auditiva. Por exemplo, para uma compo- 
nente de frequência perto dos 1000 Hz, têm sido detec- 
tadas diferenças tão notáveis como sômente 3 Hz. O 
exacto mecanismo desta assombrosa descriminação de 
frequências (a qual pareca violar a relação de incerteza 
Af.4t> 1) é ainda desconhecido. 

A maioria dos métodos de análise e sintese da voz 
são baseados na descrição paramétrica do «short-ceps- 
trum» ou de representações equivalentes, das quais é de 
salientar o logaritmo do módulo do espectro de potên- 
cia do sinal, visto que se o «short-cepstrum» permite 
conhecer com exactidão o «pitch», aquele dá-nos a 


(4) Durante o nosso trabalho temos apresentado o sinal falado filtrado entre O e 2500 Hz, o que comparado com o que aca- 
bamos de mencionar acerca da banda telefónica (200-3 200 Hz) justifica um esclarecimento, 

Na realidade, o facto de usarmos a banda O — 2500 Hz foi imposto pelo material que tinhamos à nossa disposição (frequência 
de amostragem do computador limitando o máximo da banda). Contudo o sinal falado filtrado entre O e 2500 Hz não difere sensivel- 
mente do filtrado entre 200 e 3200 Hz, apresentando igualmente inteligibilidade e boa qualdade. Especialmente no que respeita ao 
«pitch», a frequência que lhe corresponde está sempre abaixo dos 2500 Hz e no que respeita às formantes, como veremos adiante, 


os seus valores encontram-se, geralmente, abaixo dos 2500 Hz. 
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conhecer as frequências e amplitudes das formantes. 

O mais antigo método de análise-síntese da voz 
empregando uma descrição paramétrica do «short-ceps- 
trum o do «short-time» logaritmo do módulo do espectro 
de potência é «spectrum-channel vocoder». 

Neste, a envolvente espectral, espectro das for- 
mantes) H(w,t) é representado geralmente por 10 a 
20 amostras espaçadas ao longo do eixo da frequência 

O espectro da estrutura fina (sinal «avoiced») 
| P(c, t) | por um parâmetro adicional que mede o «pit- 
ch» ou a correspondente frequência fundamental Tia 
dos sons «voiced» e é igual a zero para os sons «unvoi- 
ced» ou silêncio. 

Um diagrama-bloco de um «spectrum-channel voco- 
der» é mostrado na figura 8.1. 

O sinal falado é separado em 14 bandas espectrais 
continuas com larguras de banda entre 100 e 400 Hz 
cobrindo as frequências de 200 a 3200 Hz. Esta largura 
de banda, típica em sinais telefónicos permite alta inte- 
ligibilidade e boa qualidade como já foi dito. 

A saída de cada filtro é ligada a um rectificador 
seguido de um filtro passa-baixo, à saída do qual temos 
o valor médio no tempo da amplitude do sinal resposta 
ao impulso unitário do tracto vocal para cada banda de 
frequências. 

Repare-se que o detector «voiced» — «unvoiced» 
utiliza geralmente o resultado do detector do «pitch» 
(por razões de economia de meios) ou seja, a técnica 
de «cepstrum», podendo, por exemplo, o bloco detector 
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«voiced» — «unvoiced» funcionar do seguinte modo: 
mede as amplitudes de todos os picos que aparecem no 
«cepstrum» e analisa (por comparação) se existe algum 
muito maior que todos os restantes, decidindo que o 
sinal é «voiced» em caso afirmativo e que é «unvoiced» 
em caso negativo. 

Os 14 canais referidos representam no seu conjunto 
a envolvente do «short-time» logaritmo do módulo do 
espectro do sinal falado, ou sejam as formantes. 

Um detector «voiceds — «unvoiced» e um detector 
do «pitch» fixam a estrutura fina do sinal falado e produ- 
zem um sinal correspondente de banda estreita. 

Estes 15 sinais (14 canais mais o «pitch») combi- 
nam-se num único sinal de 15 X 20 Hz. 

Como se vê, a banda de transmissão é somente 
um décimo da requerida inicialmente para transmitir 
o sinal falado. 

Chama-se a atenção que na prática é necessária 
uma banda adicional para a separação do sinal transmi- 
tido em canais individuais de sinal. 

Para uma transmissão digital, a economia ainda é 
maior. Assim para uma amostragem de 40 por segundo, 
transmitindo três «bits» por amostra para o sinal dos 
14 canais e 6 «bits» por amostra para o «pitch», a saída 
do analizador «vocoder» pode ser representada por um 
total de 40.(14 X 3 + 6) = 1920 «bits» por segundo, 
isto é, menor que um vigésimo do número de «bits» por 
segundo necessário para representar o sinal falado origi- 
nal (49 000 «bits» por segundo). 
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Contudo a qualidade da voz dos «vocoders» exis- 
tentes, funcionando com tão baixo valor de «bits» é noto- 
riamente inferior. 

No sintetisador, os sinais dos 14 canais são usiili- 
zados para modular um sinal (com espectro de frequên- 
cias plano) de frequência corresnondente ao «pitch» para 
os sons «voiced» ou ruído branco para os sons «unvoi- 
cedy», sendo em seguida filtrados pela correspondente 
banda, 

Um «vocoder» de canal esnectral como o que des- 
crevemos foi construído primeiramente por H. W. Dudley. 

Apesar da grande economia alcançada, na transmis- 
são da envolvente espectral em termos de amostras es- 
pectrais, esta ainda pode ser melhorada. A razão reside 
em que amostras espectrais, especialmente de sinais 
de canais adjacentes do «vocoder» estão significativa- 
mente correlacionadas. Assim Kramer e Mathews suge- 
riram uma transformação linear dos sinais de canal num 
pequeno número de sinais não correlacionados antes da 
transmissão [3]. Reduções de cerca de dois para um 
no número de canais têm sido realizadas na prática por 
este método. 

Um outro método para reduzir o número de canais 
activados é conhecido por «peak-picking». 

Resumidamente, consiste em transmitir os sinais 
das canais que num dado instante são maiores (tipica- 
mente um terço do seu número são máximos relativos 
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num dado instante). Assim vemos que uma redução de 
três para um pode ser conseguida. 

Para uma transmissão digital, mesmo tendo em 
conta a banda adicional necessária para transmitir os 
números dos canais seleccionados a economia é grande. 

Seguidamente, iremos indicar outros tipos de «vo- 
coders», assim como, resumidamente, as bases em que 
se fundamentam. 

Começamos pelos «vocoders» de formantes («for- 
mant-vocoders»). As envolventes espectrais de muitos 
dos sons falados são caracterizadas por máximos que re- 
presentam as ressonâncias do tracto vocal, ou seja, as 
formantes. Ora, sabemos também que abaixo dos 3000 
Hz, existem tipicamente três formantes para a voz dos 
adultos. Assim num «vocoder» de formantes o analizador 
procura detectar as frequências das formantes principais, 
assim como as suas amplitudes as quais são utilizadas 
para controlar, respectivamente, filtros variáveis e modu- 
ladores do sintetisador. A excitação do sintetisador é 
feita do mesmo modo que no «spectrum channel voco- 
der». Um esquema exemplificativo do «vocoder» de for- 
mantes com sintetisador paralelo é mostrado na fi- 
gura 8.2. 

Um modelo ainda mais económico, no que respeita 
à transmissão das formantes é o chamado «Pattern-Ma- 
tching Vocoder». Num dado instante, é determinada a 
melhor semelhança entre o «short-time» spectrum e um 
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número fixo de modelos espectrais armazenados, no que 
diz respeito às formantes. Um código representando a 
melhor semelhança é transmitido para o sintetisador, 
activando o espectro seleccionado, o qual irá modular 
o sinal excitador. 

Esquematicamente temos representado um «Pat- 
terr-Metching Vocoder» na figura 8.3. 

Filtros são incluídos para fazer uma transição suave 
de um som modelo para outro. 

Temos ainda o «vocoder» de correlação que se 


MANIA IM Ah! A Lull! ini ALiAM 
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50 


baseia, como o nome indica, na determinação da auto- 
correlação do sinal falado e os «vocoders» de Alta Fideli- 
dade de que é exemplo o «Voice-excided vocoder». 

Este último tem como função transmitir sinais fa- 
lados de alta qualidade com uma banda de 10 kHz sobre 
uma banda de 3 kHz. 

Do que se disse atrás em relação ao ouvido percebe- 
-se facilmente que o grande problema nos «vocoders» se 
põe no sintetisador com a função de excitação, ou seja, 
com o correcto conhecimento do «pitch». Por isso, no 
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Fig. 9.1 — Sinal da fig. 3.3 imerso em muído branco. Relação sinal-ruído: 3,4 dB 
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Fig. 9.2 — «Cepstrum» do sinal da fig. 9. Pico: 7 = 7,0 mseg. (Pico da origem multiplisado por zero) 
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Fig. 9.3 —- Sinal da fig. 3.3 imerso em ruído branco. Relação Sinal-Ruído: O dB 
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Fig. 9.4 — «Cepstrum» do sinal da fig 9.3 
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«voiced-excided vocoder», da banda de 10 kHz iniciais, 
uma banda baixa de 2 kHz é enviada directamente, tendo 
no sintetisador duas funções de excitação, ou seja, de 
fornecer o «pitch» o mais correctamente possível, como 
se pode ver na figura 8.4 

Como acabamos de ver, é fundamental a separação 
da envolvente da estrutura fina, ou seja, a separação das 
formantes do «pitch». 

Até à data foi o «cepstrum» o método que melhores 
resultados forneceu nesta divisão e na sua determinação. 


S — RESULTADOS DA APLICAÇÃO DO «CEPSTRUM» 
CONSOANTE A RELAÇÃO SINAL RUÍDO 


Vamos agora ver que para sinais imersos em ruído 
de modo que a relação sinal-ruído tenha pequenos va- 
lores a técnica de «cepstrum» deixa de fornecer resul- 
tados satisfatórios, apresentando um exemplo. 

Na figura 9.1 está representado o sinal da figura 
3.3 imerso em ruído branco, sendo a relação sinal-ruído 
3,4 dB. Na figura 9.2 está representado o respectivo 
«cepstrum» e como se vê para esta relação sinal-ruído 
ainda é bem nítido o pico para & = 7,0 mseg, corres- 
pondente ao «pitch», como se viu no parágrafo 3. 

Na figura 9.3, por sua vez, está representado o 
mesmo sinal imerso em ruído branco, mas desta vez a 
relação sinal-ruído é de O dB. Para este caso ainda está 
presente o pico correspondente a & = 7,0 mseg, mas 
a sua amplitude é comparável à de outros picos presentes 
(por exemplo, em torno de & = 20 e & = 40 mseg). Sem 
o conhecimento que tínhamos «a priori» de qual o valor 
«pitch», não poderíamos decidir qual dos picos corres- 
ponderia a esse «pitch». 

Para relações sinal-ruído ainda menores, o pico 
pode estar totalmente camuflado nas irregularidades do 
«cepstrum». 


10 — REFERÊNCIA A OUTROS MÉTODOS DE DETEC- 
ÇÃO DO «PITCH» E FORMANTES. 


Dado que o «cepstrum» não fornece resultados sa- 
tisfatórios para os casos de sinais imersos em ruído 
com uma pequena relação sinal-ruído e ainda para um 
certo número de outros casos como, por exemplo, os re- 
feridos no parágrafo 3, foram desenvolvidos outros mé- 
todos de detecção do «pitch», que não estudaremos, 
porque o nosso cbjectivo era o de estudar as possibili- 
dades da técnica de «cepstrum», mas aque desejamos 
fazer referência para um conhecimento mais exacto 
dessas mesmas possibilidades e para uma possível con- 
tinuação dos nossos etudos sobre os sinais falados. 


M. Noll desenvolveu três métodos que permitem 
obter resultados para aqueles casos que referimos e para 
os quais o «cepstrum» não satisfaz. São eles a determi- 
nação dc Espectro Produto das Harmónicas, do Espectro 
Soma das Harmónicas e a estimação da Máxima Veroesi- 
milhança estudados no seu trabalho que referimos com 
c número [4]. 
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Um programa para o traçado de curvas 


de nível em pequenos calculadores 


SUMÁRIO 


Indica-se um algoritmo para o traçado automático 
de isolinhas de funções cujos valores são conhecidos 
nos nós de uma malha rectangular. Embora já existam 
bastantes trabalhos sobre este problema implementados 
em grandes instalações, o autor desconhece outras ver- 
sões que sejam facilmente aplicáveis a pequenos cal- 
culadores. 

Apresenta-se a listagem do programa em BASIC 
para um calculador Hewlett-Packard 9830. Existe no 
NEEM (IST) uma versão para o Hewlett-Packard 98204A- 
-001. 


INTRODUÇÃO 


Este programa surgiu da necessidade de traçar 
automaticamente as curvas de nível de uma função cujos 
valores fossem conhecidos nos nós de uma malha rec- 
tangular. 

Foi escrito numa linguagem específica do Calcula- 
dor 9820A da Hewlett-Packard e supõe instalada a 
opção 001 (429 registos de memória) e a existência 
de um «9862A Calculator Plotter Hewlett-Packard» aco- 
plado ao Calculador. 

Posteriormente foi feita uma tradução do programa 
para BASIC pelo aluno bolseiro José Francisco Oliveira 
Santos. 

O programa foi realizado e afinado no Núcleo de 
Estudos de Engenharia Mecânica. 


DESCRIÇÃO DO ALGORITMO 
O método utilizado consiste em definir no domínio 


da função a estudar um reticulado, dando o valor da 
função nos nós da malha estabelecida. 
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J. J. AMARANTE DOS SANTOS 


ABSTRACT 


An algorithm for contour drawing of functions 
whose values are known at grid points of a rectangular 
mesh is presented. There are many implementations for 
problems of this kind for great computers but no versions 
for mini-calculators are known to the author. 

A program listing written in BASIC for a Hewlett- 
Packard 9830 Calculator is presented. Another version 
for the Hewlett-Packard 98204-001 is available at the 
NEEM (IST). 


O programa permite traçar curvas de nível que se 
fecham dentro do domínio e curvas de nível que come- 
çam numa fronteira e acabam noutra (ou nessa) fron- 
teira, 

Os nós dividir-se-ão em duas categorias: 


1) Nós em que a função tem um valor supe- 
rior ao da curva (+); 
2) Nós em que a função tem um valor inferior 
ao da curva (—). 
As malhas dividir-se-ão também em dois tipos: 


Malha simples —- Por dentro só passa 


n) o -0 
1 curva de nivel 
e----e 
+ + 
7») to---0" Malha dupla — Duas curvas de nível 
abortas ou uma fechada 
Õ---s 
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CRITÉRIOS 


O traçado das curvas de nível processa-se fazendo 
interpolações lineares entre nós consecutivos, Usam-se 
interpolações lineares porque se pensa que se o n.º de 
pontos for insuficiente, ao fazermos interpolações mais 
sofisticadas estamos a introduzir informação que nos 
possibilitará o traçado duma curva que pode nada ter 
a ver com a função em estudo em virtude da informação 
introduzida ser em grande parte arbitrária. 

O facto ds num determinado caso aparecerem 
muitas malhas do tipo 2) quer precisamente dizer que 
não se possui suficiente informação sobre a função para 
se tratar a curva com segurança. À título de exemplo 
apresentam-se para uma mesma distribuição de nós as 
curvas que se obtêm utilizando o critério a), o critério b), 
ou o critério c). 


critério b) critério c) 


critério a) 


É bem patente neste exemplo que a informação 
que se dispõe sobre a função é insuficiente e que qual- 
quer dos critérios corresponde a uma hipótese tão plau- 
sível como as cutras, sendo talvez o critério a) que 
apresenta maior arbitrarisdade, pois entre duas curvas 
de nível, em que uma está contida na outra, apresenta 
pontos em que a função toma valores acima e abaixo do 
valor quo tem na curva. Os critérios b) e c) são perfeita- 
mente equivalentes. 

Como norma geral, poder-se-á portanto dizer que 
quando a percentagem de malhas de tipo 2) for elevada, 
o gráfico não tem grande significado, tornando-se neces- 
sário um refinamento da malha. 

Para que o programa não ficasse muito complicado 
era necessário no entanto escolher um critério, que foi 
o b). Não é muito difícil alterar o programa para este 
seguir o critério c). 

Quando as curvas correspondentes ao nivel dese- 
jado estão bem definidas mas passam ambas por uma 
mesma malha, o critério b) dá uma só curva caso se tra- 
tem de duas depressões, enquanto que se forem duas 
elevações dará duas curvas. Com o critério c) passa-se 
o contrário. 


DESCRIÇÃO DO PROGRAMA 
1.1. — GENERALIDADES 


Quando o programa se inicia os valores da função 
nos nós da malha devem já estar metidos nos registos 
correspondentes. 

Por cada segmento que une dois nós consecutivos 
não pode passar mais do que uma curva de nível visto 
a interpolação ser linear. 
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Vê-se neste exemplo que a única zona mal definida 
é a que corresponde à malha do tipo 2). 


A cada segmento que une dois nós consecutivos 
da malha corresponde um determinado bit que se pode 
encontrar no estado de «aceso» ou de «apagado». 

Quando o programa se inícia, todos os bits se 
devem encontrar no estado de «apagado». 

À medida que o programa vai correndo, quando 
uma curva de nível intersecta um determinado segmento, 
o bit que corresvonde a esse segmento é posto no es- 
tado de «aceso». Este processo destina-se a evitar que 
uma mesma curva de nível seja traçada duas ou mais 
vezes. 

A lógica do programa pressupõe que os pontos da 
curva de nível a traçar se encontram sempre sobre seg- 
mentos unindo nós consecutivos. Numa terceira versão 
do programa terá que ser previsto um ciclo inicial em 
que se houver algum nó em que o valor da função seja 
igual ao da cota da curva se adicione àquele uma quan- 
tidade insignificante para o traçado da curva mas que 
evite possíveis descontroles (por exemplo 10 -Jx valor 
da função no nó). 


1.2. — VARRIMENTO DAS FRONTEIRAS 


Partindo portanto do pressuposto que os pontos 
da curva assentam sempre sobre segmentos, o programa 
começa por verificar se existe algum ponto duma curva 
du nível de cota dada no início do programa, no 1.º 
segmento da 1.º linha da malha. Se não existir, passa 
ao 2.º segmento e assim sucessivamente até esgotar a 
1.º linha. Se nesta fronteira não houver nenhum ponto 
o programa vai pesquisar nas outras fronteiras do do- 
mínio. 
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Quando é descoberto um ponto da curva numa fron- 
teira é aceso o bit correspondente a esse segmento e 
inicia-se o traçado duma curva de nível cuja evolução 
se passa a descrever. 
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Supondo que o ponto detectado se situava no seg- 
mento (2-3) o programa iria em 1.º lugar testar o seg- 
mento (3-7) para ver se passava por lá alguma curva 
de nível; Em caso afirmativo uniria o ponto a com o 
ponto do segmento (3-7), seguindo assim o critério b) 
no caso da malha ser dupla. Caso não existisse nenhum 
ponto com a cota dada no segmento (3-7) iria pesquisar 
no segmento (2-6) e depois em (6-7). A ordem de pes- 
quisa de (2-6) e (6-7) é irrelevante para o traçado da 
curva: O que realmente interessa para a execução do 
critério b) é que se 3 for um ponto (+) o 1.º segmento 
e pesquisar seja (3-7). 

Dada a sequência como a pesquisa é feita mal o 
programa detecte um ponto num daqueles 3 segmentos, 
une-o logo com o ponto a. Se por exemplo existisse um 


o: 
Entrada da cota da curva 


cíclo que apa- 


ga os bits 


varrimento da 
fronteira 


nto fronteíro a seguir ao ínicíal 


varrimento 
dos pontos * 
interiores 


Recomeça os testes no segme 
no segmento interior a seguir ao inicial 


46: 


Chegou ao se nto 
' GRRRE 


Chegou a uma fronteira 


bit aceso, 


42: 


Aceso 


ponto em (3-7), depons do traçado do 1. segmento de 
recta da curva de nível o programa iria pesquisar a exis- 
tência de pontos em (3-4), (7-8) e finalmente (4-8), re- 
petindo-se este processo até que a curva de nível inter- 
sectasso um outro segmento de recta pertencendo a uma 
fronteira. 

Ao longo do traçado da curva o programa vai acen- 
dendo os bits correspondentes aos «segmentos de recta 
unindo nós consecutivos» que a curva vai intersectando. 


Quando a curva chega a uma fronteira a pena é 
levantada e recomeça a pesquisa de pontos ao longo das 
fronteiras no segmento imediatamente a seguir âquele 
em que tinha começado o traçado da curva. Será opor- 
tuno acrescentar, que em virtude do programa testar o 
estado do bit correspondente a um segmento em que se 
detectou um ponto da curva, quando o programa passar 
pelo segmento em que terminou a 1.º curva de nível tra- 
çada, o traçado desta já não será repetido. 

Nota: em virtude do segmento em que começa 
uma curva de nível já não voltar a ser testado, é irrele- 
vante acender o bit que lhe corresponde. 


1.3. — VARRIMENTO DOS PONTOS INTERIORES 


Para detectar curvas de nível fechadas são var- 
ridos os segmentos interiores Só é necessário varrer um 


Exíste (valor da cota compreendido entre os valores 


dos extremos do segmento) 


19; 


apagado 


fronteí 


Existe 


Apagado 


22: 


Inicialização 


(Guarda-se o segmento 
inicial 


Traçado da curva 
(descrito na pág. 5) 


Outro caso (vai pf 47:) 


tipo de segmentos (horizontais ou verticais). 30 DISP «NUM. DE PONTOS SEGUNDO OX, OY»; 


Quando num segmento interior (segmento que não 40 INPUT X1,X2 
pertence a uma fronteira) é detectado um ponto duma 50 GOSUB 2270 
curva começa-se o traçado da mesma maneira que para 60 R9 = SQOR( (X1 + 1)*X2 — 2) 
as curvas começando numa fronteira. 70 IF R9 = INT (R9) THEN 90 
O traçado da curva pára quando se passa pela 2, 80 R9 = INT (R9) + 1 
vez pelo segmento inicial. 90 R8= 12 + 2*ºR9 
Como o processo de «acender bits» não permite 100 A1 = FNI (0) 
que se «acenda» o mesmo bit duas vezes, quando se 110 DISP «COTAS: 


inicia o traçado de uma curva fechada não se acende 
o bit do segmento inicial, vindo aquele a ser aceso 
quando se fecha a curva. 
Depois do varrimento dos pontos interiores é pedido 
outro valor para a cota. 150 IF RO F [XY] THEN 170 
Nota: O programa não permite que haja mais do 160 FIX, Y] = 1.0001ºF [X, Y] 
que 1296 nós (o que corresponde a 36 X 36) mas isso 170 NEXT X 


120 INPUT RO 
130 FOR Y = 1 TO X2 
140 FOR X=1 TO X1 


não constitui verdadeiramente uma limitação pois pode- 180 NEXT Y 
-se sempre dividir o domínio da função em vários sub- 190 A = 12 | 
domínios. 200 R [12 + R9] = 1 


210 A=A+1 
220 R[A] = O 


TESTE EFECTUADO 230 R[A+R9] = 2*R [A+ R9 — 1] 
240 IF A <= 11 + R9 THEN 210 
Na figura abaixo reproduzem-se as curvas obtidas 250 I8 = 0 
para uma função f (x, y) que era nula em todos os nós 260 Y = 1 
excepto nos definidos pelas coordenadas: 270 X = 1 


(2,1), (3,1); (2,3); (3,4); (5,3), (5,4), (6,3), (6,4); 
(1,5), (1,6); (3,6), (4,6), (4,7); (8,y = 1,11) 
em que tinha o valor 5, 


280 B=X+1-I8 

290 C=Y+I8 

300 GOSUB 1960 

310 IF 15 =1 THEN 460 

20 IF (Y=1) AND (X=X1 — 1) THEN 390 

30 IF (Y = X2) AND (X=X1 —1) THEN 410 
340 IF (X=1) AND (Y=X2-1) THEN 430 

FIGURA DU a50 IF (X=X1) AND (Y=X2-1) THEN 630 


sd Es 360 X = B 
370 Y=C 
380 GOTO 280 


ww 


390 Y = X2 
É> 400 GOTO 270 
410 8 = 1 
420 GOTO 260 
| ==" 430 X =X1 
MtiiaeemSS 440 = 1 
450 GOTO 280 
As curvas traçadas correspondem às cotas 4, 3, 460 7 =1 


261. 
O tempo para traçar cada uma das curvas foi de 
aproximadamente 1 minuto e 10 segundos. Este tempo 


470 GOSUB 2020 
480 IF |7 = 1 THEN 320 


: 490 R4 =X 
que parece um pouco excessivo deve-se em grande parte 
à existência de muitas curvas em virtude do exemplo 500 R6 = B 
ser um teste. Numa função melhor definida do que esta 910 R5 =Y 
existirá para cada cota muito menos curvas. 920 R7=C 
Esse facto reduzirá o tempo pois para o programa 530 IF 18 = 1 THEN 590 
se aperceber da inexistência de curvas de cota 6 demo- 540 13 = 0 
rou apenas 24 segundos. 550 IF (Y = 1) AND (I8= 0) THEN 570 


560 IF ((Y = X2) AND (I8 = 0)) OR ((X = 1) AND 
(18 = 1) ) THEN 610 


LISTAGEM DO PROGRAMA 570 14 = 0 
580 GOTO 830 

10 DIM R [250], FS [36,36] 590 13 = 1 

20 19 = 0 600 GOTO 550 
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Synopsis of articles published in «Técnica» n.º 428 


C. D. U. 621,3:517.512.2 
IRIARTE JOSÉ ARAUJO ESTEVES 
e 
JOSE ALBERTO TEIXEIRA 


«CEPSTRUM» APLICATION TO SPEECH ANALYSIS 
«VOCODERS» 


This paper, reprinted from a technical report of Laboratório 
do Análise de Sonais, |. À. C., starts with a short description 
of speech sjgnals and a brief theoretical analvsis of ceps- 
trum. These techniques are then applied to the determination 
of speech parameters. The papsr ende with an introduction 
to vocoders. 


J. J. AMARANTE DOS SANTOS 


A PROGRAM FOR CONTOUR DRAWING 
IN MINI-CALCULATORS 


An algorithm for contour drawing of functjons whose values 
aro known at grid points of a rectangular mesh is presented. 
There are many implementations for problems of this kind for 
great computers but no versjons for mini-calculators are 
known to the author. 

& program listing written in BASIC for a Hewlett-Packard 
9830 Calculator is presented. Another version for the Hew- 
let-Packard 98204-001 is available at the NEEM (IST). 
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U. D. C. 517.(084):681.14 


J. LAGINHA SERAFIM 
e 
MARCIANO L. GUERREIRO 


TESTS AND OBSERVED BEHAVIOUR OF ROCK 
FOUNDATIONS OF DAMS 


The paper presents results of many vyear's efforts to interpret 
deformabjlity tests of foundation rocks of concrete dams. 
Results of «in situ» deformability tests depend on the direc- 
tion, depth of tests and thegeological details of the site. Cor- 
relatjons with depth and dip and strike of the formations are 
very clear. Methods of interpretation used and comparisons 
with laboratory results are discussed and results of many 
tests at site in various rocks are presented. New findings 
on the effects of grouting in deformability are indicated and 
tho most recent information is given on rocks foundatjon 
behaviour, based on observation with long length strain me- 
ters (extensometers). 


U. D. C. 628.5:551.5104 
J. J. DELGADO DOMINGOS 
M. TERESA PONTES 
AN AIR POLLUTION GAUSSIAN DIFFUSION MODEL 


This report presents a Gaussian dispersion model for gassous 
pollutants which accounts for the effects of multiple ponctual 
sources, temperatura inversion, altmospheric stability and 
variatjons in both wind speed and direction. 

The stability of the atmosphere is classified into Pasquill's 
classes and the siandard deviatjons of the concentration 
lav. oz) are calculated by means of a numerical approxi- 
matior to Pasquill-Gitford's plots. 

A few examples ilustrating the effect of atmospheric stabjlity 
and oí inversion on ground level concentration and of the 
interference of nearby sources are also worked out. 
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C. D. U. 624.131.388:627.82 
J. LAGINHA SERAFIM 
E] 
MARCIANO L. GUERREIRO 


ENSAIOS DE FUNDAÇÕES ROCHOSAS DE BARRAGENS 
E SEU COMPORTAMENTO 


Este trabalho apresenta os resultados dos esforços de vários 
anos para interpretar ensajos de deformabilidade de funda- 
ções rochosas de barragens de betão. Os resultados de en- 
saios de deformabilidade «in situ» dependem da direcção e 
profundidade dos pontos ensaiados e das condições geoló- 
gicas, sendo muito claras as correlações com a profundidade 
e inclinação das formações. São discutidos os métodos de 
interpretação usados e feitas comparações com resultados 
obtidos em laboratório, sendo anresentados os resultados de 
diversos ensaios «in situ» em diferentes rochas. Indicam-se 
novos factos sobre os efeitos da jniecção na deformidade 
agora encontrados e dão-se informações recentes sobre o 
comportamento das fundações em rocha, baseadas em obser- 
vações com medidores da extensões de grande comprimento 
(extensometros) 


C. D. U. 628.5:551.510.4 
SJ. J. DELGADO DOMINGOS 
E 
M. TERESA PONTES 


UM MODELO GAUSSIANO PARA A DISPERSÃO 
DE POLUENTES NA ATMOSFERA 


Apresenta-se um modelo Gaussiano de dispersão de poluen- 
tes na atmosfera, considerando o efeito de fontes múltiplas 
pontuais, a existência de inversão de temperatura em alti- 
tude, a variabilidada da direcção e velocidade do vento e da 
estabilidade da atmosfera. 

A estabilidade da atmosfera é classificada segundo o mé- 
todo da Pasquill sendo o desvio padrão da concentração de 
poluentes calculado a partir de aproximações numéricas de 
gráficos de Pasquil-Gifford. 

No final são apresentados alguns exemplos que ilustram o 
efeito da variação da estabilidade da atmosfera, da existência 
de inversão e a interferência de fontes próximas. 


Ano L — Julho 1975 
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C. D. U. 621.3:517.512.2 


IRIARTE JOSE ARAUJO ESTEVES 
a 
JOSÉ ALBERTO TEIXEIRA 


«CEPSTRUM» 
ESTUDO EXPERIMENTAL DA SUA APLICAÇÃO AO SINAL 
FALADO, «VOCODERS» 


Esto trabalho, que constitui relatório de um estágio realizado 
no Laboratório de Análise de Sinais do |. À. €., é iniciado 
com uma análise resumida do sinal falado e um breve estudo 
teórico da técnica do cepstrum. Segujdamente é apresentada 
a aplicação desta técnica na determinação de alguns parã- 
metros de sinais daquele tipo. Por úúltimo, faz-se uma jntro- 
dução ao estudo dos vocoders. 


C. D. U. 517.5/084):681,14 
J. J. AMARANTE DOS SANTOS 


UM PROCRAMA PARA O TRAÇADO DE CURVAS 
DE NIVEL EM PEQUENOS CALCULADORES 


Indica-so um algoritmo vara o traçado automático de isoli- 
nhas de funções cuios valores são conhecidos nos nós de 
uma malha rectangular. Embora iã existam bastantes traba- 
lhos sobro este problema implementados em grandes insta- 
lações, o autor desconhece outras versões que seiam facil- 
mento aplicáveis a pequenos calculadores. 

Apresenta-se a listagem do oroorama em BASIC para um 
calculador Hewlett-Packard 9830, Existe no NEEM (IST) 
uma versão para o Hewlett-Packard 98204-001. 


= 0000 =D = 


2... ,.—————. — — 
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1800 
1810 
1820 
1830 
-FLX, 
1840 
1850 
1860 
1870 
1880 
1890 
1900 
1910 
1920 
1930 
1940 
1950 
1960 
1970 


1980 
1990 
2000 
2010 
2020 
2030 
2040 
2050 
2060 
2070 
2080 
2090 
2100 
2110 
2120 
2130 
2140 


414 


4 = 1 

RETURN 

END 

R3= (RO — F [X,Y]) /(FI[B,C] 
Y]) 

IF I9 = O THEN 1900 
IF X=B THEN 1930 
B=X+R3* (B—X) 
PLOT BC 

I9 = O 

RETURN 

7 =0 

GOSUB 2020 

GOTO 1850 
C=>7YEFRPIC-T; 
GOTO 1870 

END 

I5 = 0 


IF (F [X,Y) — RO)J* (F [BC] — RO) <= 0 


THEN 1990 

RETURN 

5=1 

RETURN 

END 

IF B/X THEN 2200 

IF (X521) AND (X&X1) THEN 2250 
SO =(X1-1)"X2+INT((CA+Y)/2) 
IF X4/X1 THEN 2070 
S0=S0+X2-1 

S1 = INT (SO/R9) 

SO = SO — S1*R9 

IF SO = O THEN 2220 

S1=13+51 

IF I7=1 THEN 2140 

R [SI] =R[S1] +R [12 + R9 + SO] 
RETURN 

82 = À [81] 


2150 
2160 
2170 
2180 
2190 
2200 
2210 
2220 
2230 
2240 
2250 
2260 
2270 
2280 
2290 
2300 
2310 
2320 
2330 
2340 
2350 
2360 
2370 
2380 
2390 
2400 
2410 
2420 
2430 
[5,3] 
2440 
[4,6] 
2450 
2460 
2470 
2480 
2490 


S2 = 52/R [12 + R9 + 50] 
Z = INT (S2)/2 

IF Z-INT (Z) THEN 2250 
7=0 

RETURN 


SO =(X1 = 1 *[Y — 1) INT L(X + B)/2) 


GOTO 2070 
So = R9 

81=81-1 

GOTO 2100 

RETURN 

END 

SCALE 1,X1,1,%2 
DISP «QUER EIXOS? SIM> 
INPUT L1 

IF L1 =0 THEN 2350 
XAXIS 1,1 

YAXIS 1,1 

XAXIS X2,1 

YAXIS X1,1 

RETURN 

END 

DEF FNI (Z) 

FOR X=1 TO 11 
FOR Y=1 TO 11 
FIX,Y1=0 

NEXT Y 

NEXT X 


: NÃO >0»; 


FIiZAI = FS =F[2,3) = FIA) =E 


= F[54]=F[6,3]=5 


F[6,4] = F[1,5] = F[1,6] = F[3,6] =F 


= F[4,7]=5 

FOR Y=1 TO 11 
F [8 y]=5 

NEXT Y 

RETURN O 

END 
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610 
620 
630 
640 
650 
660 
670 
680 
690 
700 
710 
720 
730 
740 
750 
760 
770 
780 
790 
800 
810 
820 
830 
540 
850 
860 
870 
880 
890 
900 
910 
920 
930 
940 
950 
960 
970 
980 
990 
1000 
1010 
1020 
1030 
1040 
1050 
1060 
1070 
1080 
1090 
1100 
1110 
1120 
1130 
1140 
1150 


1160 


1170 


1180 


TEC 


4 =1 

GOTO 830 
I8 = O 

= ee 
X=1 
B=X+1 
C=Y 
GOSUB 1960 


IF [5 = 1 THEN 780 
IF X o X1-1 THEN 760 
N ==4 
IF Y <= X2-2 THEN 740 
GOTO 110 
Y=Yy+1 
GOTO 660 
X=5 
GOTO 660 
= 
GOSUB 2020 
IF |I7 = 1 THEN 700 
6.=4 
GOTO 490 
GOSUB 1830 
A=(-1) 14 
IF |3=1 THEN 1010 
C=Y+AÃ 
IF F[X,Y]>FI[X+1,Y] THEN 930 
GOSUB 1300 
IF 15 =1 THEN 1150 
GOSUB 1450 
IF 15 =1 THEN 150 
GOTO 970 
GOSUB 1450 
IF [5 =1 THEN 1150 
GOSUB 1300 
IF 15 =1 THEN 1150 
Y= c 
B=X+41 
GOSUB 1830 
GOTO 1150 
B=X-a 
IF FIX,YI>EFIX,Y+ 1] THEN 1080 
GOSUB 1570 
IF |5=1 THEN 1150 
GOSUB 1710 
IF 15 = 1 THEN 1150 
GOTO 1120 E 
GOSUB 1710 
IF |5=1 THEN 1150 
GOSUB 570 
IF 15 =1 THEN 1150 
C=Y+1 
X =:B 
GOSUB 1830 
IF (X = R4) AND (Y = R5) AND (I3 = |8) THEN 
1190 


IF ((X= 1) OR (X=X1)) AND (13 = 1) THEN 
1230 

IF ((Y=1) OR (Y =X2)) AND (I3=0) THEN 
1230 

GOTO 840 
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1190 
1200 
1210 
1220 
1230 
1240 
1250 
1260 
1270 
1280 
1290 
1300 
1310 
1320 
1330 
1340 
1350 
1360 
1370 
1380 
1390 
1400 
1410 
1420 
1430 
1440 
1450 
1460 


B-—- R6 
C=-R7 
PEN 

GOTO 700 
B— R6 
C=R7 

XxX — R4 

Y =R5 
PEN 

GOTO 320 
END 
X=X+1 
B-—- X 
GOSUB 1960 


IF |5 = 1 THEN 1366 
IF |5 = 1 THEN 1380 


GOTO 1420 
GOSUB 1830 
GOTO 1340 
Y=Y-I4 
3 =1 

4 =1 
RETURN 
X=X—1 
RETURN 

END 

B=X 
GOSUB 1960 


1470 IF |5 = 1 THEN 150€ 


1480 
1490 
1500 
1510 
1520 
1530 
1540 
1550 
1560 
1570 
1580 
1590 
1600 
1610 
1620 
1630 
1640 
1650 
1660 
1670 
1680 
1690 
1700 
1710 
1720 
1730 
1740 
1750 
1760 
1770 
1780 
1790 


IF 15 = 1 THEN 1529 
RETURN 

GOSUB 1830 

GOTO 1480 
Y=Y-—I4 

I3=1 

4 = O 

RETURN 

END 

TENTA 

C=Y 

GOSUB 960 

IF I5 = 1 THEN 164€ 
IF I5 = 1 THEN 1660 
EV =4 


RETURN 
GOSUB 1830 
GOTO 1610 
A=K-— +14 
I3=0 

I4 = O 
RETURN 

END 

C=Y 

GOSUB 1960 
IF I5 = 1 THEN 1760 


IF I5 = 1 THEN 1780 
RETURN 

GOSUB 1830 

GOTO 1740 
XK=X—1+14 
I3=0 
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